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ABSTRACT 

A technique for the use of a digital computer for simul- 
taneous generation of feedback control system transient 
response and sensitivity to parameter variations is developed. 
The sensitivity model for a basic linear control system is a 
cascaded replica of the original system with the sensitivity 
pickoff points corresponding to the error and feedback signals 
of the original system. This is convenient for computer anal- 
ysis with DSL/360. For systems with nonlinear elements, the 
sensitivity model is implemented by making the instantaneous 
gain equal to the slope of the nonlinear function. Sensi- 
tivity analysis is utilized for improvement of system 
performance, determining parameter tolerances, and predicting 
changes in the response due to parameter variations. 
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INTRODUCTION 



I . 



The early development of sensitivity analysis was 
applied to circuit theory where it is useful to know how 
sensitive the response of the circuit is to variations of 
circuit parameters such as resistance, capacitance, and 
inductance. There are many different definitions of sensi- 
tivity available. H. Bode [1] defines a logarithmic or 
normalized sensitivity function of a variable, x^(t, a j ) , 



Hortowitz [2] adapted the inverse of this definition which 

is widely used. J. Leeds [3] presented sensitivity function 

for incremental changes in parameters and S. R. Parker [4] 

recently generalized the technique of sensitivity analysis for 

nonlinear circuits. Much of this work entailed the definition 

of R. Tomovic [5] who defined the sensitivity coefficients as 
x. 3x. x, 3x. 

i _ i 

S = -r . Kokotovic [6] adapted the definition S = tt-s 

a^ 3aj ^ a, 3in 

for use in control systems and developed a method of generat- 
ing sensitivity coefficients using analog computers. In this 
thesis the method of Kokotovic [61 is expanded for use with 
digital computer techniques utilizing block diagrams and 
transfer functions of forward and feedback loops in control 
systems. The IBM/360 Digital Simulation Language (DSL) 
program [7] is very well suited for modeling a control system 
and obtaining the time response and sensitivity coefficients 
simultaneous ly . 





and I . M. 



a 



j 



3 Jin x^ 
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DSL is a continuous system simulator whose non-procedural, 
problem-oriented input language has the attractive building- 
block approach of digital-analog simulators with the added 
power of logical and algebraic equation notation. The DSL 
system provides a basic set of 35 blocks from which a phys- 
ical system may be modeled. Included in this set are such 
conventional analog components as integrators, comparators, 
and transfer functions. The Kokotovic [6] method of sensi- 
tivity analysis is readily adaptable to the DSL program as 
control systems in block-diagram form can be programed in 
the s -domain. 

Section II of this thesis applies the Kokotovic [6] 
method of generating sensitivity coefficients of the basic 
feedback control system for the forward and feedback gain 
parameters. Sensitivity functions are derived and applied 
to a third-order system. Examples are given for the use 
of sensitivity coefficients in analysis, parameter adjust- 
ment for improvement of system performance, and parameter 
tolerance considerations. Section III applies the 
Kokotovic [6] method of sensitivity analysis to multiloop 
systems. Gain and time-constant sensitivity coefficients 
are derived and applied to a regulated power generator as an 
example. Section IV demonstrates the Kokotovic [6] method 
for the general case of differential equations, and Section 
V is devoted to the nonlinear system. 

Sensitivity analysis of a control system offers an 
analytical method for determining those parameters of a 



10 



dynamic system which have the greatest influence on the 
system performance. High-speed digital computers give 
very precise solutions for the system equations, but it is 
impossible to realize exact mathematical models and para- 
meter values. Additional variations in parameters may also 
be caused by wear, aging, shifting of operating point, 
temperature, and etc. It is beneficial to know the effect 
of these small parameter variations on system performance. 
Sensitivity analysis can give this necessary insight so 
that all parameters of a dynamic system are not given equal 
importance. This is more efficient than repetitive computer 
solutions with parameters changed by increments. Analysis 
utilizing sensitivity coefficients can put emphasis on 
specific component tolerances and points that are critical 
to the performance specifications. 
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II. BASIC CONTROL SYSTEM SENSITIVITY 



A. INTRODUCTION 

The Kokotovic [6] method for computing sensitivity 
coefficients for various parameters is restricted to linear, 
time-invariant, deterministic systems. Although it is 
true that physical systems do not have these ideal properties, 
linear approximations for these systems have proven adequate 
in many useful cases. When parameter values are to be 
approximated, it is important to estimate the influence of 
each parameter on the system response. This can be done 
accurately by sensitivity analysis. Also, limits of indi- 
vidual parameter variations may be deduced and tolerance 
levels formulated. 

The Kokotovic [6] method of solving for sensitivity 
coefficients has the advantage that all of the functions 
are generated simultaneously with the time-domain solution. 
Kokotovic [6] gives analytic proof that this is possible. 

As an engineering approach to analysis, it is common to 
approximate a high-order system by an equivalent second- 
order system to enable using graphs and tables for second- 
order response characteristics. Qualitative observations 
of graphs of sensitivity coefficients can give insight into 
the transient performance in terms of overshoot, settling 
time, rise time, and etc., for higher degree characteristic 
equations as will be seen in examples. 
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B. GAIN CONSTANT SENSITIVITY COEFFICIENTS 



The sensitivity coefficient S„ is defined as the ratio 

K 

of the change in output to the fractional change in the 
parameter K. 



_ 3C 3C 

K " 3K/K 3£nK 



( 2 . 1 ) 



where C is the controlled output and K is the parameter being 
investigated. 

The canonical form of a feedback control system after 
Laplace transformation is shown in Fig. 2.1. The definitions 
of terms are as follows: 



K 

KG (s) 

F 

FH (s) 
KG ( s) FH (s) 
R(s) 
E 
B 



direct gain constant 
direct transfer function 
feedback gain constant 
feedback transfer function 
loop transfer function 
reference input 
error signal 
feedback signal 



The basic equations are as follows: 



C, > = KG (s) 

R ; 1+KG ( s ) FH ( s ) 



E, . _ 1 

R' S 1+KG ( s ) FH ( s ) 



B, = KG (s) FH (s) 

R V ' 1+KG ( s ) FH ( s ) 



(2.2a) 

(2.2b) 



(2.2c) 



In the work that follows, capital letters denote the 
Laplace transformation of the functions described except for 
the letters K and F which are constants . 
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R + 
>- 



^ E V 


KG ( s ) 


C 


1 












B 


FH ( s ) 











Fig. 2.1. Canonical form of a feedback, control system. 




Fig. 2. 2. Feedback control system and sensitivity model. 
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1.5 




TIME (SEC) 



Fig. 2. 3. Basic control system response and sensitivity 
for K = 50. 
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Applying (2.1) to (2.2a) yields the gain sensitivity 
coefficients 



RKG 



C 



S 



K 



( 1+KGFH) 2 



1+KGFH 



(2.3) 



S 



-RK 2 G 2 FH 



C (-KGFH) 
1+KGFH 



(2.4) 



F 



(1+KGFH) 2 



These two gain sensitivity functions are shown in block 
diagram in Fig. 2.2. This shows that the sensitivity model 
is an exact replica of the original system with the controlled 
output as the input to the sensitivity model. The error and 
feedback signals' of the sensitivity model are the desired 
sensitivity functions. It should be noted that the sensi- 
tivity model requires no new components or transfer functions 
and all quantities are generated simultaneously. 

C. APPLICATION 

In this thesis the equations for the feedback systems 
are programed on the IBM 360 computer using the Digital 
Simulation Language (DSL) program for both the transient 
solution and the sensitivity coefficients. This program 
has the advantages of programing in the s~domain in trans- 
fer function form and the accuracy of the digital computer 
with graphical and numerical results. A sample program is 
included in Appendix A and the numerical solution and 
sensitivity coefficients are listed in Table 2.1. Fig. 2.3 
is the graph of the output and sensitivity coefficients. 
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The direct transfer function is 



KG 



50 

s ( s+4) (s+6 ) 



(2.5) 



For example, data from Table 2.1 when S K is a maximum 

Time T = 1.0 

Gain K = 50 

Sensitivity S = 0.70863 

K 

Input R = unit-step 

Output C = 0.96456 

For large variations (2.1) becomes 

SK = ^ (2.6) 

Let the gain, K, change 10%, then 

ATC 

AC - S R = 0.70863(0.1) = 0.070863 



This shows that the output will change approximately + 0.071 

if the gain varies from 50 to 55 or 45. This can be 

verified from Table 2.2 for the following data: 

K = 55 

T = 1.0 

C = 1.0331 

C' = C + AC = 0.96456 + 0.070863 = 1.03541 

and for 

K = 45 

T = 1.0 

C = 0.89133 

C' = C-AC = 0.96456 - 0.07086 = 0.89370 
This shows that sensitivity coefficients for changes of gain 
as large as 10% produce accuracy within 0.3% when compared to 
incremental intervals as denoted in partial derivatives. 
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Another example for the use of sensitivity coefficients 
is to compute the gain required to reduce the peak overshoot 
of the system transient response to a more desirable specifi- 
cation of one-half its initial value. Table 2.1 gives the 
following data: 



K = 50 

Peak overshoot = 0.2 

T = 1.5 

S R = 0.36568 

Rearranging (2.6) 



AK 



ACK 

S K 



K' 



K - AK 
50 (1 



= K(l-g^) 
b K 

0.1 

0.36568 



) 



36 



(2.7) 

( 2 . 8 ) 



Table 2.2 for K = 36 shows that the peak overshoot is reduced 
to 0.1 and Fig. 2.4 shows the reduction of peak overshoot by 
a factor of one-half. 

Any other design specification, such as settling time or 
rise time, can be adjusted by utilization of the sensitivity 
coefficient. In the example it should be noted that reducing 
the overshoot by reduction in gain, increases the rise time 
by 0 . 4 second. 

Tolerance considerations for parameters are encompassed 
within the definition of sensitivity coefficient. 



Rearranging (2.6) 

Ar 

Tolerance = AK/K = — 

b K 



(2.9) 
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1.2 



1.0 



0.8 



0.6 



0 .4 



0.2 



C.O 




O.C 2.0 4.0 6.0 8.0 10.0 

TIME ( SEC ) 



Fig. 2. 4. Overshoot reduction of 0.1 by utilization of 
sensitivity coefficients. 
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0.0 2.0 4.0 6.0 8.0 10. 

TIME ( SEC ) 



Fig-2. 5- Gain tolerance of 14$ for overshoot of ±0.05. 
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BASIC SERVO SENSITIVITY COEFFICIENTS 



TIME 


OUT 


SENSK 


0.0 




0.0 


0.0 


1 .OOOOE- 


■01 


6.518 2 E-03 


6.516CE-C3 


2.00C0E- 


■01 


4. 1385E-02 


4. 1258E-02 


3 .0000 E- 


■01 


1.1 1 65E-0 1 


1. 1054E-C1 


4.0300E- 


■01 


2 • 1292 E-0 1 


2 .08 12 E-C 1 


5.0000E- 


•01 


3. 3623E-01 


3. 2212E-C 1 


6 .0000 E- 


■01 


4 . 71 47 E-0 1 


4 . 3885 E-C l 


7.0000E- 


•01 


6 . 09 1 1 E-0 1 


5.4526 F-C 1 


8.00C0E- 


■01 


7.41C5E-C1 


6. 3042E-01 


9 .0000 E- 


■01 


8 . 6 100 E-0 1 


6.864CE-C 1 


l.OOOOE 


00 


9.6456E-01 


7. 0863E-C 1 


1.1000E 


00 


1 .049 1 E 00 


6. 959CE-0 1 


1.2000E 


00 


1.1137E 00 


6. 5005 E-C 1 


1.3000E 


00 


1.1586E 00 


5. 75 42 E-0 1 


1.40C0E 


00 


1.1853E 00 


4.7816E-01 


1 .5000E 


00 


1.1959E 00 


3. 6558E-01 


1.6000E 


00 


1 . 193 1 E 00 


2. 4539E-01 


1.7000E 


00 


1.1797E 00 


1 .25 1C E-0 1 


1.8000E 


00 


1.1587E 00 


1. 1496E-C 2 


1 .9000E 


00 


1.1328E 00 


-8. 9767E-02 


2.0000E 


00 


1.1046E 00 


- 1. 7441 E-0 1 


2. 1000E 


00 


1.0760E 00 


-2. 3964E-01 


2.2000E 


00 


1.0489E 00 


-2. 84 18 F- 01 


2.3OC0E 


00 


1.0245E 00 


-3.0811 E-C 1 


2.4000E 


00 


1.0036E 00 


-3. 1273E-C 1 


2.5000E 


00 


9. 8672 E-0 1 


- 3. 00 32 E-0 1 


2.6000E 


00 


9. 7404E-0 1 


-2. 7383E-01 


2.7000E 


00 


9.6544E-01 


-2. 3663E-C1 


2.8O00E 


00 


9 . 6056 E-0 1 


- 1 . 9225 E-C 1 


2.900CE 


00 


9. 5895E-9 1 


-1 . 44 14 E-C 1 


3.0000E 


00 


9 . 600 2 E-0 1 


-9. 5469E-C2 


3. 10C0E 


00 


9.631 8 E-0 1 


-4. 8996 F-C 2 


3.2000F 


00 


9.6782E-01 


-6.9429E-03 


3.300CE 


00 


9.7336E-01 


2 .90 47 E-0 2 


3.4000E 


00 


9. 7932E-0 1 


5. 79 14 E-C 2 


3.5000E 


00 


9. 8525E-0 1 


7. 9165 E-C 2 


3.600C-E 


00 


9.9C83E-01 


9.280 8E-C 2 


3.7000E 


00 


9.9581 E-0 1 


9. 9281E-C2 


3.8000E 


00 


l.OOOOE 00 


9.9353F-C2 


3.9C0CE 


00 


1.0C34E 00 


9. 40 37 E-C 2 


4.0000E 


00 


1.0058E 00 


8.4486F-C2 


4.1000E 


00 


1.0075E OC 


7 . 1906 E-C 2 


4.2000E 


00 


1.0084E 00 


5. 7483E-02 


4.3000E 


00 


1.0086E 00 


4. 2310E-C2 


4.4000E 


00 


1.0083E 00 


2.735C E-C2 


4.5000E 


00 


1.0075E 00 


1.33 90 E-0 2 


4.6000E 


00 


1.0065E 00 


1.0386E-03 


4.7000E 


00 


1.0053E 00 


-9. 2916E-03 


4.8000E 


00 


1.0041E 00 


-1. 7368E-02 


4.9000E 


00 


1 .0028 E 00 


-2.3124E-n2 


5.0000E 


00 


1.0017E 00 


-2. 6629E-02 


5.100CE 


00 


1.0007E 00 


-2.8C75F-02 


5.2000E 


00 


9 • 9984E-0 1 


-2 . 7730 E-C 2 


5.3000E 


00 


9. 9918E-01 


-2. 5919F-C2 


5.4000E 


00 


9 • 9869 E-0 1 


-2 . 2996E-02 


Table 2.1. Computer output data for unity 



■f 



SENSF 

- 0 . 0 

-2. 2129E-06 
-1.2685E-C4 
-1. 1099E-C3 
-4.7931E-03 
- 1 • 4109E-02 
-3. 2619E-02 
—6 • 3846E— 02 
-1. 1063E-0 1 
-1. 7460E-01 
-2.5593E-01 
-3. 5322E-01 
-4. 6363E-0 1 
-5. 8319E-01 
-7. 0712E-C1 
-8.3030E-01 
-9. 4767E-01 
-1.0546E CO 
-1. 1472F 00 
-1.2226E 00 
-1.2790E 00 
-1.3157E 00 
-1.333 IE OC 
-1.3326E 00 
-1.3163F 00 
-1.2870E OC 
-1.2479E 00 
-1.2021E 00 
-1. 1528E 00 
-1 . 1031 E 00 
-1.0555E OC 
-1.0122E 00 
-9. 7476E-01 
-9.4432E-01 
-9 . 2140 E-C 1 
-9. 0609E-01 
-8 • 9802E-0 1 
-8. 9653E-01 
-9.0O66E-01 
-9.0933E-C1 
-9. 2136E-01 
-9. 3558E-01 
-9. 5088E-01 
-9. 6627E-C1 
-9. 8091E-01 
-9.9414E-C1 
-1.0055E 00 
-1.0146E OC 
-1.0214E 00 
-1.0260E OC 
-1.0283E 00 
-1.0288E ^0 
-1.0276E 00 
-1.0251E nn 
-1.0217F 00 



feedback ay stem. 
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BASIC SERVO SENSITIVITY COEFFICIENTS 



Table 2 



TIME 

5.5^00E 00 
5.6000E 00 
5.7COOE 00 
5.8000E 00 
5.9000E 00 
6.0000E 00 
6.1000E 00 
6.200CE 00 
6.3000E 00 
6.4000E 00 
6.5C00E 00 
6.6000E 00 
6.7000E CO 
6.8000E 00 
6.9C00E 00 
7.0000E 00 
7.I000E 00 
7.200CE 00 
7.3000E 00 
7.4900E OO 
7.500CE 00 
7.6000E 00 
7.7C00E 00 
7.8000E 00 
7.9C0CE 00 
8.0000E 00 
8.100CE 00 
8.2000E 00 
8.3C00E 00 
8.40COE 00 
8.5000E 00 
8.600CE 00 
8.70C0E 00 
8.8000E 00 
8.9000E 00 
9.0000E 00 
9.I000E 00 
9.2000E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000E 00 
9.7C00E 00 
9 .8000 E 00 
9.900CE 00 
l.OOOOE 01 



OUT 

9.9838E-01 
9.9823 E-0 1 
9.9821E-01 
9.9829E-01 
9.9846E-01 
9. 9868E-0 1 
9. 9893E-0 1 
9.9920E-01 
9 . 9945E-0 1 
9.9969E-01 
9.9989 E-0 1 
1 .000 1 E 00 
1.0002E 00 
1.0003E 00 
1.0003E 00 
1.0C04E 00 
1.0004E 00 
1.0003E 00 
1.0003E 00 
1.0003E 00 
1 .0002 E 00 
1.0002E 00 
1 .000 I E 00 
1 .0001 E 00 
l.OOOOE 00 
9.9998 E-0 1 
9 . 9995E-0 1 
9.9994E-01 
9 . 9992 E-0 1 
9.9992E-01 
9.9992 E-0 1 
9.9993E-01 
9. 99 9 3 E-0 1 
9. 9994 E-0 1 
9.9995E-01 
9. 9997E-0 1 
9. 9998 E-0 1 
9 . 9999E-0 1 
9.9999E-0 1 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OCOOE 00 
l.OOOOE 00 



-1.9314E-02 
-1. 52C8E-02 
-1.098CE-02 
-6. 8853E-03 
-3.1281 E-03 
1 • 40 61 E-04 
2.82 57E-03 
4. 8800 E-03 
6. 3002E-C3 
7.11 94E-C3 
7.3986E-03 
7.2 1 71 E-03 
6 • 6686E-03 
5. 8481 E-03 
4. 8504E-03 
3.7631 F-03 
2 • 6627E-C 3 
1.6125E-C 3 
6. 6310 E-04 
-1. 5068E-04 
-8.08 72 E-04 
-1 . 3018E-C3 
-1.6327E-03 
-1.8120 E-03 
-1. 8587E-03 
-1.7906 E-03 
-1.6360 E-03 
-1.4186E-03 
-1. 1616E-C3 
-8. 8680E-C4 
-6. 1244E-C4 
-3.5483E-C4 
- 1 . 2463E-C 4 
6. 9797 E-0 5 
2.2489E-C4 
3. 3873E-C4 
4. 1258E-C4 
4.4990E-C4 
4 • 5526E-04 
4. 3404 E-04 
3.923 8E-C4 
3.3724E-04 
2.7323E-04 
2.0635E-04 
1. 40 67 E-04 
8. G466F-05 



SENSF 

-1.0177E 00 
-1.0134E 00 
-1.0092E 00 
-1.0052E 00 
-1.0016E 00 
-9. 9854E-C 1 
-9.9611E-01 
-9 . 9432 E-0 1 
-9.931 5E-01 
-9. 9257E-01 
-9. 9249E-01 
-9. 92 8 4 E- 01 
-9. 9352E-0 1 
-9. 9444E-01 
-9.9549E-01 
-9. 9661 E-0 1 
-9. 9771E-01 
-9. 9873E-C 1 
-9 . 9965 E-0 1 
-1.0004E 00 
-1.0010E 00 
-1.0015E 00 
-1.0017E 00 
-1.0019E CO 
-1.0019E OC 
-1.0C18E OC 
-1.0016E Of 
-1.0014E 00 
-1.001 IE 00 
-1.0008E 00 
-1.0005E 00 
-1.0C03E OC 

-l.ooo i e no 

-9. 9987E-01 
-9. 9973E-C1 
-9. 9963E-0 1 
-9. 9956E-0 1 
-9. 9954 E-0 1 
-9. 9954E-C 1 
-9 • 9957E-0 1 
-9. 996 1 E— 0 1 
-9.9967E-01 
-9. 9974E-0 1 
-9. 9981E-01 
-9 . 9987 E-0 1 
-9. 9993E-01 



. 1 . ( Continued ) 
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TIMF 

0.0 

1 .0000 E-C 1 

2 . COOOE-O 1 

3. C000E-01 
4.0000E-01 
5. OOOOE-O 1 
6.0000E-01 
7 • OOOOE-O 1 
P. 00 00 E-01 
9.0CC0E-01 
l.OOOOE CO 
1.1C00E 00 
1.2000E 00 
1.3000E 00 
1.4000E 00 
1.5000E 00 
1.60C0E 00 
1.700CE 00 
1.800CE 00 
1.9C00E 00 
2.0000E 00 
2.1000E 00 
2.2000E 00 
2.3000E 0" 
2.40C0E 00 
2.50C0E 00 
2.6C00E 00 
2.7000E 00 
2.8000E 00 
2.9000E 00 
3.0000E 00 
3.1000E 00 
3.2000E 00 
3.3000E 00 
3.40C0E 00 
3.5000E 00 
3 • 6000 E CO 
3.70C0E 00 
3 .80 DOE 00 
3.9000E OC 
4.0000E 00 
4.1000E 00 
4.2000E 00 
4.3000E 00 
4.4000E 00 
4.5000E 00 
4.6000E 00 
4.7000E 00 
4.80C0E 00 
4.9000E 00 
5.0000E 00 
5.1CCCE 00 
5.2000E 00 
5.3000E 00 
5.4000E 00 



Table 2.2. 



K =50 K=55 K=45 



K=3 6 



C.O 

6.5 182E-03 
4. 1385E-02 
1. 1165E-01 
2 . 1292E-01 
3 • 3623E-01 
4.7147E-01 
6.091 IE-01 
7.4105E-31 
8.6100E-01 
9. 6456E-01 
1.0491E 30 
1.1137E 00 
1.1586E 00 
1.1853E 00 
1.1959E 00 
1 . 193 1 E 00 
1.1797E 00 
1.1587E 30 
1.1328E 03 
1.1046E 30 
1.0760E 30 
1.0489E 00 
1.0245E 00 
1.0036E 00 
9. 8672E-0 1 
9.7404E-01 
9. 6544 E- 3 1 
9.6056E-01 
9.5895E-01 
9. 6002E-0 1 
9.6318 E-0 1 
9. 6782E-01 
9. 7336E-0 1 
9.793 2 E-0 1 
9. B525E-01 
9.9083E-31 
9 . 95 8 1 E-0 1 
l.OOOOE 00 
1.0034E 00 
1.0058E 30 
1.0C75E 00 
1.0084E 00 
1.0C86E 00 
1 .0083 E 30 
1.0075E 30 
1.0065E 00 
1.0053E 00 
1.0G41E 00 
1.0028E 30 
1.0017E 00 
1.0007E 30 
9.9984E-01 
9.99 18E-01 
9 . 9869E-0 1 



OUT 

0.0 

7 • 1698 E-0 3 
4.5510E-02 
1.2269E-31 
2.3369E-31 
3. 6830E-01 
5. 1503E-01 
6.6301 E-01 
8.0304E-91 
9 . 2801 E-0 1 
1 .0331 E 30 
1.1156E 00 
1.1747E 30 
1.2114E 00 
1.2276E 00 
1.2264E 00 
1.2112E 30 
1.1857E 00 
1.1536E 00 
1.1180E 00 
1.0823E OC 
1.0478E 00 
1.0173E 00 
9.9169E-01 
9.71 57E-3 1 
9.5715 E-0 1 
9.4823E-01 
9 . 4430 E- 3 1 
9 • 4463E-0 1 
9.4837E-01 
9. 5462E-0 1 
9. 6253E-31 
9.7123E-31 
9. 8001 E-01 
9. 8836E-01 
9.958 1E-0 1 
1.0021E 30 
1.0C70E 00 
1.0105E 30 
1.0127E 00 
1.0136E 30 
1.0135E 00 
1.0126E 30 
1.0111E 30 
1.0C92E 00 
1.007GE 30 
1.0C49E 30 
1.0028E 00 
1.0013E 30 
9. 994 8 E-01 
9. 9828E-01 
9 . 9742E-3 1 
9 . 9689 E-0 1 
9. 9666E-01 
9.9669E-01 



OUT 

0.0 

5. 8666E-03 
3.7258E-02 
1.0053E-01 
1 .9206E-01 
3.03 87 E-0 1 
4.2726E-31 
5 . 5396E-0 1 
6 . 7694E-0 1 
7.9071E-01 
8.9133E-01 
9. 7634E-01 
1.0446E 00 
1.0961E 00 
1.1317E 00 
1.1529E 30 
1.1617E 00 
1 . 160 IE 30 
1.1505E 30 
1.1350E 00 
1.1158E 00 
1.0946E 00 
1.0729E 30 
1.0520E 00 
1.0328E 00 
1.0159E 00 
1 .00 18 E 00 
9.9052E-01 
9.8215E-01 
9 * 7648 E-0 1 
9. 7324E-01 
9.72C8E-01 
9.7262E-01 
9. 7449E-01 
9 . 7729 E-0 1 
9. 8063E-0 1 
9.8437E-01 
9.8810E-01 
9. 9166E-01 
9. 9490 E-01 
9. 9773E-0 1 
1.0001E 00 
1.0019E 90 
1.0033E 30 
1.0042E 00 
1.0047E 30 
1.0048E 00 
1.0047E 00 
1.0043E 00 
1.0038E 00 
1 .0032 E 00 
1.0026E 00 
1.0019E 00 
1.0013E 00 
1.0008E 00 



OUT 

0.0 

4. 6935E-03 
2.9823E-02 
8 .06 10E-02 
1 • 5427E-0 1 
2.4494E-01 
3 . 46C7E-9 1 
4. 5154E-01 
5 • 5613E-01 
6.5572E-01 
7.4727E-01 
8.2873E-01 
8.9895E-01 
9 • 5749E-0 1 
1.0C45E 00 
1.0407E 00 
1.0668E OC 
1.0842E 00 
1.0941E 30 
1.0977E 00 
1.0963E 30 
1 .091 IE 00 
1.0832E 00 
1.0736E 00 
1.0630E 30 
1.0522E 00 
1.0416E 00 
1.0316E 00 
1.0226E OC 
1.0146E 00 
1.0C78E 09 
1.0023E 30 
9.9790E-01 
9 . 9459 E-0 1 
9. 9227E-01 
9.9081 E-01 
9. 9007E-91 
9 • 8992E-0 1 
9.9023E-01 
9.9089E-01 
9. 9179E-0 1 
9. 9283E-01 
9 . 9394E-01 
9.9506 E-0 1 
9.9613E-01 
9.9712E-01 
9.9801 E-0 1 
9. 9878E-0 1 
9 • 9943E-0 1 
9. 9996E-01 
1.0004E 30 
1.0007E 30 
1.0C09E 00 
1.0010E 00 
1.0010E 00 



Computer output data for K = 55, 5°, ^5, and J>6. 
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5.5o5ie E oo 
5.6000E 00 
5.7000E 00 
5.8000E 00 
5.9000E 00 
6.0C0CE 00 
6.1000E 00 
6.2000E 00 
6.3000E 00 
6.4000E 00 
6.5000E 00 
6.6000E 00 
6.700CE 00 
6.8000E 00 
6.9000E 00 
7.0000E 00 
7.1000E 00 
7.2000E 00 
7.300CE 00 
7.40Q0E 00 
7.5000E 00 
7.60C0E 00 
7.7000E 00 
7.8000E 00 
7.9000E 00 
8.0000E 00 
8.1000E 00 
8.2000E 00 
8.3000E 00 
8.4000E 00 
8.5000E 00 
8.6000E 00 
8.7000E 00 
8.8C00E 00 
8.9000E 00 
9.0000E 00 
9.1000E 00 
9.20C0E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000 E 00 
9.7000E 00 
9.8000E 00 
9.9000E 00 
1.0000E 01 



K=50 

9.983(Se-01 
9.9823E-01 
9.9821E-01 
9. 9829E-01 
9.9846E-01 
9. 9868E-0 1 
9.9893E-01 
9.9920 E-0 1 
9.9945E-01 
9.9969E-01 
9.9989E-01 
1.0001E 00 
1.0002E 00 
1.0003E 00 
1.0003E 00 
1.0004E 00 
1.0004E 00 
1.0003E 00 
1.0003E 00 
1.0003E 00 
1 .0002 E 00 
1.0002E 00 
1.0001E 00 
1.0001E 00 
l.OGOOE 00 
9 • 9998 E-0 1 

9. 9995 E-0 1 
9.9994E-01 
9.9992E-01 
9.9992E-01 
9.9992E-01 
9.9993E-01 
9. 9993E-0 1 
9.9994E-01 
9.9995E-01 

9.9997 E-0 1 

9 .9998 E-0 1 
9.9999E-01 
9.9999E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 



K=55 

9. 96§¥e-0 1 
9.9729E-01 
9.9776E-01 
9.9828E-01 
9.9881 E-0 1 
9. 993 1 E-0 1 
9.9975E-0 1 
1.0001E 00 
1.0004E 00 
1 .0006 E 00 
1.0008E 00 
1.0008E 00 
1.0008E 00 
1.0008E 00 
1.0007E 00 
1.0005E 00 
1.0004E 00 
1 .0003 E 00 
1.0002E 00 
1 .000 IE 00 
9 • 9997 E-0 1 
9. 9989E-0 1 
9 • 9984E-0 1 
9.9901 E-0 1 
9. 9980E-01 
9 • 9980 E-0 1 
9.9981E-01 
9 . 9984E-0 1 
9 • 9986 E-0 1 
9. 9989E-01 
9.9993E-01 

9. 9996 E-0 1 
9. 9998E-0 1 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 90 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 



K=45 

OUT 

1.0003E 00 
9. 9992E-0 1 
9.9962E-01 
9.9940E-01 
9 . 9926E-Q 1 
9.9919E-01 
9.9917E-91 
9.9920 E-0 1 
9.9927E-01 
9.9936E-91 
9 .9946 E~0 1 
9.9957E-01 
9.9968E-01 
9.9978E-01 
9. 9988E-01 
9.9996E-01 
l.OOOOE 00 
1.0001E 00 
1.0001E 00 
1.0C01E 00 
1.0001E 00 
1.0001E 00 
1.0001E 00 
1 .000 IE 00 
1.0001E 00 
1.0001E 00 
1.0C01E 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
9 .9999 E-0 1 
9.9998E-01 
9. 9998E-0 1 
9.9997E-01 

9.9997 E-0 1 
9 • 9997E-0 1 

9.9997 E-0 1 
9. 9998E-01 

9.9998 E-0 1 

9. 9998 E-0 1 
9.9999E-01 
9.9999E-01 
9. 9999E-01 
9.9999E-01 
l.OOOOE 00 



Table 2.2. (Continued) 
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OUT 

1.0010E 00 
1.0010E 00 
1.0009E 00 
1.0008E 00 
1.0007E 00 
1.0006E 00 
1.0005E 00 
1.0004E 00 
1.0003E 00 
1.0002E 00 
1.0001E 00 
l.OOOOE 00 
9.9998E-01 
9.9995E-01 
9 • 9992 E-0 1 
9. 9990E-01 

9.9989 E-0 1 
9. 9989E-01 
9. 9989E-01 

9 . 9990 E-0 1 
9.999 1 E-0 1 
9.9992E-01 
9. 9993E-0 1 
9.9994E-01 
9.9995E-0 1 
9. 9996 E-0 1 
9. 9997E-0 1 

9 .9998 E-0 1 
9.9999E-01 
9.9999E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 

i.ooooe on 
l.OOOOE 00 
1 .OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
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If the design criteria is to restrict the peak overshoot to 
+ 0.05 then (2.9) gives 



Figure 2.5 shows the results of the extreme values of K = 43 
and K. -= 57 which give ’the tolerable • limits of overshoot. 



A. INTRODUCTION 

The Kokotovic [ 6 ] method of sensitivity analysis is 
expanded to investigate gain and time-contant variations in 
multiloop systems using a’ procedure similar to the one 
developed for the basic feedback control system. This proce- 
dure is developed using,-the block diagram of a regulated 
power generator which incl-udes an amplifier-exciter, a series 
RC compensating circuit in the feedback loop, and a load as 
shown in Fig. 3.1. 

B. GAIN-CONSTANT SENSITIVITY COEFFICIENTS 

Terms are defined as follows: 

KiGi(s) 5 Direct transfer function for the amplifier- 



AK/K 



0.05 
0 . 366 



14% 



III. MULTILOOP CONTROL SYSTEM SENSITIVITY 



exciter * 



K 2 G , 2 ( s ) = Transfer function for the generator 



K 3 G 3 ( s ) = Transfer function for the load 



FH (s) 



Feedback v. ,.mt 



or transfer function 



The basic equations are as follows: 




KG. ( s ) K 2 G , (s)K 3 G's (s) 



1+KxG, is) EH ( s) ^KjGi (s) K 2 G 2 (s) K 3 G 3 (s) 



(3.1a) 
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e 2 

r" 



(s) 



1 

1+K ! G i (s)FH(s)+KiGi (s)K 2 G 2 (s)K 3 G 3 (s) 



(3.1b) 



Ei 

R- (s > 



1+K i G i (s) FH (s) 

1+K i G i ( s ) FH ( s) +K i G ! (s)K 2 G 2 (s)K 3 G 3 (s) 



(3.1c) 



B i 

I" (S) 



KiGi (s)FH(s) 

1+K i G i (s) FH (s) +KiGi (s) K 2 G 2 (s) K 3 G 3 (s) 



(3. Id) 



As before, capital letters denote the Laplace transform 
of the function described, except for K and F which are gain 
constants, and T which is a time constant. 

To obtain the system sensitivity function apply (2.1) 
to (3.1a) to give 



S 



K, 



RKiGiK 2 G 2 K 3 G 3 

(l+KiGiFH+KiGxKzGzKaGs) 2 

(3.2) 



C 

1+KiGiFH+KiGiK 2 G 2 K 3 G 3 

This sensitivity function corresponds to the error signal, 

E 2 , in the original multiloop system but with the controlled 
output, C, as the input to the system. This can be seen by 
comparing (3.2) to (3.1b) with C replacing R. As in the 
basic control system of Section II, the sensitivity model for 
gain parameters can be realized by cascading two identical 
systems and using the error signals in the sensitivity model 
as the desired gain sensitivity coefficients as shown in 
Fig. 3.2. 



The load-gain sensitivity function is derived by applying 
(2.1) to (3.1a), yielding 
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Fig. 3.1. Multiloop system. 



RK1G1K2G2K3G3 ( 1 +KiGiFH) 



S K 3 = (I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 ) 

C (I+K 1 G 1 FH) 
I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 



(3.3) 



This sensitivity function corresponds to the error signal, 
Ej , as can be seen by comparing (3.3) to (3.1c) with C 
replacing R as the input. 

For the feedback-gain sensitivity function apply (2.1) 
to (3.1a), yielding 



RK1G1K2G2K3G3 (-K1G1FH) 
I+K1G1FH+K1G1K2G2K3G3) 2 

C (-KiGjFH) 

1 +K 1 G 1 FH+K 1G1K2G2K3G3 



(3.4) 



This sensitivity function corresponds to the negative of the 
feedback signal, Bi, as can be seen by comparing (3.4) to 
(3. Id) with C replacing R as before. 

By applying (2.1) to the system transfer function, 
sensitivity functions have been derived for the various gain 
parameters and the pickoff points for the sensitivity 
functions can be selected by inspection as shown in Fig. 3.2. 
No new components have been required and all sensitivity 
functions can be generated simultaneously. 
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Fig. 3. 2. Multiloop sensitivity model. 



c. 



TIME-CONSTANT SENSITIVITY COEFFICIENTS 



Since time -cons tan t variations can also affect the 
response, sensitivity coefficients for LR and CR network 
combinations can add helpful information to the analysis 
of a modeled system. For this derivation, let 

C KG _ KG (Ts+1) , 

R 1+KGFH 1+Ts (KGF+1) 



where H = 
circuit. 



= — —r and T is the time constant of an RC feedback 
Ts + 1 

Applying (2.1) to (3.5) yields 
-R (KG) 2 FTs 

S = 

T [1-t-Ts (KGF+1) ] 2 



(3.6) 



C(-KGFTs) 1 

[1+Ts (KGF+1) ] Ts+1 



and the feedback-gain sensitivity is 

„ _ C(-KGFTs) 

b F 1+Ts (KGF+1) 

Therefore, by combining (3.5) and (3.6) 
S T = S F ^Ts+1 ^ 



(3.7) 



(3.8) 



This indicates that the RC time-constant sensitivity can 
be generated by using the feedback-gain sensitivity as the 
input to a block with elementary transfer function as shown 
in Fig . 3.2. 

In the direct path there are time-constants associated 
with the generator and load inductance and resistance 
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combinations. The effect of variations of these time- 
constants will add information to the analysis of the system. 
For this derivation, let 



c = KG = K 

R 1+KGFH 1+Ts+KFH v * ’ 

where G = m . and T is the time-constant of a direct path 
Ts + 1 

LR component. Applying (2.1) to (3.9) yields 



„ _ -RKTs 

b T (1+Ts+KFH) 2 

C (-Ts) 

(1+Ts+KFH) 

and the direct-gain sensitivity is 

q = C ( Ts+1 ) 

K (1+Ts+KFH) 

Therefore, by combining (3.10) and (3.11) 



(3.10) 



(3.11) 



S T = S K <fi?T> <3 - 12) 

As with the feedback RC time-constant sensitivity, only one 
elementary block has to be added in series with the direct- 
gain sensitivity model to complete the sensitivity model as 
shown in Fig. 3.2. 

D. APPLICATION 



One of the basic problems of preliminary design analysis 
is to estimate the influence of each design parameter on the 
system response and decide how improvement of the system 
can be assured. Sensitivity coefficients can be helpful in 
solving this problem. To illustrate this procedure, a basic 
set of design parameters for a regulated generator is taken to 
be amplifier gain, exciter time-constant, load-gain, 
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1.5 




TIME ( SEC ) 



Fi6*3.3. Multiloop system response and sensitivity 
functions . 



32 



0.1 



0.05 



C.O 



-0.05 



-C.l 




TIME ( SEC ) 

Fig. 3. 4. Amplifier-exciter sensitivity functions. 




0 2 A 6 8 10 

TIME (SEC ) 



Fig.3*5* Multiloop response for = 40 and = 60. 
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TIKE(SLC) 

Fig. 3 .6. Load and feedback gain sensitivity functions. 




0 2 4 6 8 10 

TIME l SitC) 

Fig. 3*7* Multiloop response for K^ = 1.1 and F= 2.2. 
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feedback-gain, and a compensator time-constant. Figure 3.3 
is a computer graph of the time response to a unit-step 
input and the five design-parameter sensitivity coefficients. 

At a glance, it is apparent from Fig. 3.3 that Ki and 

T i have relatively little influence on the response. Table 

3.1 shows the maximum value of is only 0.05 which is less 

& i 

than 10% the value of S„ . Figure 3.4 shows and S n 

K. 3 K i 1 1 

where the scale has been increased by 10 to show the general 
shape of the curves. Figure 3.5 shows the time response for 
Kj = 40 and Kj =60. Even for this 50% increase, there is 
little change in the performance. Since the response is 
very "insensitive" to Kj and Tj , little improvement can be 
achieved by adjusting these parameters and their tolerances 
can be large in relation to the other parameters. It should 
be noted that the system is practically insensitive to small 
variations of resistance in the exciter because the time- 
constant and gain both either increase or decrease by the 
same amount. The sensitivities have the same relative 
amplitude but differ in sign and therefore the effect on the 
time response is canceled. 

In Fig. 3.6 the sensitivity coefficients for the load- 
gain and feedback-gain constants are plotted. Table 3.2 
shows that their magnitudes are almost equal but are 180 
degrees out of phase. Figure 3.7 is a graph of the time 
response with their original values and the response for a 
10% increase in both parameters. Although the response is 
sensitive to such changes, the plots are almost superimposed 
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due to the one being the negative of the other. This could 
be favorable equalization if aging or temperature cause 
these gain parameters to vary by the same percentage. 

Figure 3.3 shows that the feedback-gain and the compen- 
sator RC time-constant have almost the same sensitivity but 
the time-constant sensitivity function is delayed by one 
second. 

Since fractional parameter variations are considered, 
the permissible percentage of a parameter change can be 
determined from the allowable system response variation. 

It has been shown that any general form of a control 
system that can be transformed into a system of primary 
loops can be analyzed in the same fashion as developed in 
this section. The sensitivity to variations in direct- 
and feedback-gain constants are available as error or feed- 
back signals in the sensitivity model and time-constant 
sensitivities require addition of elementary transfer 
functions as shown in Fig. 3.2. 
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Tikis, 



0.0 

1.0000E-01 
2.0000E-01 
3.0000E-01 
4.0000E-01 
5.0000E-01 
6.0000E-01 
7.0000E-01 
8*0000E“0l 
9.0000E-01 
1.0000E 00 
1* 1000E 00 
1.2000E 00 
1.3000E 00 
1.4000E 00 
1.5000E 00 
1.6000E 00 
1 • 7000E 00 
1.8000E 00 
1.9000E 00 
2* OOOOE 00 
2.1000E 00 
2.2000E 00 
2.3000E 00 
2.4000E 00 
2.5000E 00 
2.6000E 00 
2.7000E 00 
2.8000E 00 
2.9000E 00 
3.0000E 00 
3.1000E 00 
3.2000E 00 
3.3000E 00 
3.4000E 00 
3.5000E 00 
3.6000E 00 
3.7000E 00 
3.8000E 00 
3.9000E 00 
4.0000E 00 
4.1000E 00 
4.2000E 00 
4.3000E 00 
4* 4000E 00 
4.5000E 00 
4.6000E 00 
4.7000E 00 
4.8000E 00 
4.9000E 00 
5.0000E 00 
5 • 1000E 00 
5.2000E 00 
5.3000E 00 
5.4000E 00 



s Ki 



s Ti 



0.0 

5.2083E-04 
2.9189E-03 
7.4409E-03 
1.4014E-02 
2.2527E-02 
3.2866E-02 
4.4918E-02 
5* 8570E-02 
f •3713E-02 
9.0238E-02 
W0804E-01 
1.2701E-01 
1» 4705E-01 
1 •6807E-01 
1.8995E-01 
2.1262E-01 
2.3597E-01 
2.5993E-01 
2.8440E-01 
3.0930E-01 
3.3455E-01 
3.6007E-01 
3.8580E-01 
4# 1166E-0 1 
4.3759E-01 
4.6351E-01 
4.8937E-01 
5.1512E-01 
5.4069E-0 
5.6603E 
5.9111E 
6* 1586E-01 
6.4025E-01 
6.6424E-01 
6.8779E-01 
7.1086E-01 
7*3344E-01 
7.5548E-01 
7.7696E-01 
7.9786E-01 
8.1816E-0 
8.3784E-0 
8.5688E-0 
8.7528E-0 
8.9301E-0 
9.1006E-01 
9.2644E-01 
9.4214E-01 
9.5715E-01 
9.7147E-01 
9.8510E-01 
9.9804E-01 
1.0103E 00 
1.0219E 00 



:*ui 

:“01 

:-0l 

:-0l 



l 



0*0 

3.2895E-04 
1.2793E-03 
2.4687E-03 
3.7582E-03 
5.1151E-03 
6.5283E-03 
7.9903E-03 
‘.4944E-03 
•1034E-02 
1.2603E-02 
1.4194E-02 
1.5801E-02 
l*7420E-02 
1.9043E-02 
2.0664E-02 
2.2280E-02 
2.3883E-02 
2.5470E-02 
2.7035E-02 
2.8574E-02 
3.0082E-02 
3.1556E-02 
3.2991E-02 
3.4385E-02 
3.5733E-02 
3.7032E-02 
3.8281E-02 
3.9476E-02 
4.0615E-02 
4.1696E-02 
4.2717E-02 
4.3677E-02 
4.4574E-02 
4. 5408E-02 
4.6178E-02 
4.6881E-02 
4.7520E-02 
4.8093E-02 
4.8600E-02 
4.9042E-02 
4.9419E-02 
4* 9731E-02 
4.9980E-02 
5.0168E-02 
5.0293E-02 
5*0357E-02 
5.0362E-02 
5.0309E-02 
5.0203E-02 
5.0039E-02 
4.9825E-02 
4.9561E-02 
4.9247E-02 
4*8887E“02 



0.0 

-3.2650E-04 
-1* 2578E-03 
-2.4017E-03 
-3.6161E-03 
-4. 8670E-03 
-6.1426E-03 
-7.4349E-03 
-8.7368E-03 
-1.0042E-02 
"1* 1343E-02 
-1.2635E-02 
-1.3910E-02 
-1.5165E-02 
-1.6393E-02 
-1.7590E-02 
•1. 8751E-02 
-1.9872E-02 
-2.0948E-02 
- 2 » 1977E-02 
-2.2954E-02 
-2.3877E-02 
-2.4742E-02 
-2.5549E-02 
-2* 6294E-02 
-2.6976E-02 
-2.7594E-02 
-2.8146E-02 
-2.8631E-02 
-2.9048E-02 
-2.9399E-02 
-2.9681E-02 
-2. 9896E-02 
-3.0044E-02 
-3* 0126E-02 
■3.0142E-02 
-3.0092E-02 
-2.9980E-02 
-2.9805E-02 
-2.9570E-02 
■2.9276E-02 
'2* 8925E-02 
■2* 8520E-02 
■2. 8062E-02 
-2.7553E-02 
-2.6997E-02 
-2.6393E-02 
-2. 5746E-02 
•2* 5059E-02 
•2.4335E-02 
-2.3572E-02 
■2.2779E-02 
-2.1955E-02 
■2. 1103E-02 
•2* 0226E-02 



Table 3.1a. Computer output data for C, and S^. 



3 ? 



TIME' 

5.5006E00 
5.6000E 00 
5.7000E 00 
5.8000E 00 
5.9000E 00 
6.0000E 00 
6.1000E 00 
6.2000E 00 
6.3000E 00 
6.4000E 00 
6.5000E 00 
6.6000E 00 
6.7000E 00 
6. 8000E 00 
•9000E 00 
00 
00 
00 
00 



\ 



• OOOOE 

• lOOC 



_ 0 
7.2000 
7.3000 
7.4000E 00 
7.5000E 00 
7.6000E 00 
7.7000E 00 
7.8000E 00 
7.9000E 00 
8. OOOOE 00 
8.1000c 00 
8.2000E 00 
8.3000E 00 
8.4000E 00 
8.5000E 00 
8.6000E 00 
8.7000E 00 
8.8000E 00 
8.9000E 00 
9.0000E 00 
9.1000E 00 
9.2000E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000E 00 
9.7000E 00 
9.8000E 00 
9.9000E 00 
l.OOOOE 0 
1.0100E 0 
1.0200E 0 
1.0300E 0 
1.0400E 0 
1.0500E 0 
1.0600E 0 
•0700E 0 
1.0800E 0 
1.0900E 0 



• 0328C 
•0430L 
• 0526E 00 
•0615E 00 



00 

00 




.0845E 00 
• 0910E 00 
•0969E 00 

• 1021 E 00 
•1069E 00 
. 1 11 IE 00 
.1148E 00 
•1179E 00 
•1206E 00 
• 1228E 00 
•1246E 00 
• 1259E 00 
• 1268E 00 
• 1273E 00 
•1274E 00 
•1271E 00 
. 1266E 00 
• 1256E 00 
• 1244E 00 
• 1229E 00 
• 1211E 00 
• 1 191E 00 
•1168E 00 
•1143E 00 
•1116E 00 
•1087E 00 
•1056E 00 
•1024E 00 
• 0991E 00 
•0956E 00 

• 092 OE 00 
• 0883E 00 
•0845E 00 
•0806E 00 
•0767E 00 
•0728E 00 
• 0688E 00 
•0647E 00 
•0607E 00 
. 0567E 00 
•0527E 00 
• 0486E 00 
•0446E 00 
•0407E 00 
• 0368E 00 
•0329E 00 

• 0291 1 00 
• 0253E 00 
•0216E 00 



a Ki 

4. 8484E-02 
4.8035E-02 
4, 7550E-02 
4.7024E-02 
4.6467E-02 
4. 5875E-02 
4.5251E-02 
4.4601E-02 
4. 3923E-02 
4.3223E-02 
4.2501E-02 
4.1759E-02 
4.0999E-02 
4. 0227E-02 
3.9443E-02 
3.8649E-02 
3.7846E-02 
3.7040E-02 
3.6225E-02 
3.5413E-02 
3.4599E-02 
3. 3787E-02 
3.2978E-02 
3.2178E-02 
3. 1382E-02 
3.0596E-02 
2. 9821E-02 
2.9058E-02 
2.8305E-02 
2. 7569E-02 
2.6848E-02 
2.6145E-02 
2. 5457E-02 
2.4791E-02 
2. 4144E-02 
2.3515E-02 
2.2910E-02 
2.2327E-02 
2.1766E-02 
2.1228E-02 
2.0713E-02 
2.0222E-02 
1.9754E-02 
1.9313E-02 
1.8897E-02 
1.8503E-02 
1.8134E-02 
1. 7792E-02 
1.7472E-02 
1.7179E-02 
1.6908E-02 
1.6664E-02 
I.6441E-02 
L.6244E-0f 



a Ti 

-1.9329E-02 

-1.8408E-02 

-1.7475E-02 

-h 11111-81 

-1.4597E-02 
-1.3620E-02 
-1.2641E-02 
-1. 1660E-02 
-1.0681E-02 
-t9. 7037E-03 
-8.7310E-03 
-7.7655E-03 
-6.8106E-03 
-5. 8688E-03 
-4.9392E-03 
-4* 0244E-03 
-3.1285E-03 
-2.2473E-03 
-1.3897E-03 
-5. 5103E-04 
2. 6421E-04 
1.0564E-03 
1. 8213E-03 



2. 5623E-03 
3.2745E-03 
3.9590E-03 
4.6156E-03 
5. 2448E-03 
5.8424E-03 
6.4102E-03 
6.9460E-03 
7. 4535E-03 
7. 9278E-03 
8.3706E-03 
8.7855E-03 
9. 1657E-03 
9.5148E-03 
9.8346E-03 

1.0610E-02 
1. 0810E-02 
1.0979E-02 
1.1118E-02 
1. 1233E-02 
1.1320E-02 
1.1378E-02 
1. 1413E-02 
421E-02 
407E-02 
. 366E-02 

305|-0r 
L221E-0J 



h: 



1.6069E-02 1. 1 116E-02 



Table 3.1a. (Continued) 
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-vi c 
-01 7 

31 1 



00 

00 

00 

00 

00 

00 



TIME 

0.0 

1.0000E-01 
2.0000E-01 
3.0000E-0 
4.0000E 
5.0000E 
6.0000E-01 
7. 0000E-01 
8.0000E-01 
9.0000E-01 
1.0000E 00 
1.1000E 
1.2000E 
1.3000E 
1.4000E 
1.5000E 
1.6000E 
1.7000E 00 
1.8000E 00 
1.9000E 00 
2.0000E 00 
2.1000E 00 
2.2000E 00 
2.3000E 00 
2.4000E 00 
2.5000E 00 
2.6000E 00 
2.7000E 
2.8000E 
2.9000E 
3.0000E 
3.1000E 
3.2000E 
3.3000E 
3.4000E 
3.5000E 
3.6000E 
3.7000E 
3.8000E 
3.9000E 
4.0000E 
4.1000E 00 
4.2000E 00 
4.3000E 00 
4.4000E 00 
4.5000E 00 
4.6000E 00 
4.7000E 00 
4.8000E 00 
4.9000E 00 
5.0000E 



00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 



5.1000E 

5.2000E 

5.3000E 

5.4000E 



00 

00 

00 

00 

00 



0.0 

5.2081E-04 
2.9185E-03 
7,440 8E-03 
.4014E-02 
•2527E-02 
3.2866E-02 
4.4918E-02 
5.8570E-02 
7.3713E-02 
9.0238E-02 
1.0804E-01 
1,270 1E-0 1 
1.4705E-01 
1.6807E-01 
1.8995E-01 
2.1262E-01 
2.3597E-01 
2.5993E-01 
2, 8440E-01 
3.0930E-01 
3, 3455E-01 
3.6007E-01 
3.8580E-01 
4.1166E-01 
4.3759E-01 
4.6351E-01 
4, 8937E-01 
5, 1512E-01 
5.4069E-01 
5.6603E-01 
5.9111E-01 
6.1586E-01 
6.4025E-01 
6.6424E-01 
6.8779E-01 
7.1086E-01 
7, 3344E-01 
7.5548E-01 
7.7696E-01 
7.9786E-01 
8.1816E-01 
8.3784E-01 
8.5688E-01 
8.7528E-01 
8.9301E-01 
9, 1006E-01 
9.2644E-01 
9.4214E-01 
9.571 5E-01 
9, 7147E-01 
9.8510E-01 
9.9804E-01 
1.0103E 00 
1.0219E 00 



SK3 

- 0,0 

5.2079E-04 
2.9177E-03 
7, 4349E-03 

3. 27 18E-02 
4.4627E-02 
5. 8053E-02 
7.2863E-02 
8.8919E-02 
1.0609E-01 
1.2424E-0 
1.4323E-0 
1. 6293E-0 
1.8323E-01 
2.0398E-0 
2, 2507E-0 
2.4638E-0 
2.6780E-01 
2.8921E-01 
3, 1050E-01 
3,31 58E-01 
3.5235E-01 
3.7271E-01 
3.9257E-01 
4.1186E-01 
4.3051E-01 
4.4842E-01 
4.6555E-01 
4.8183E-01 
4.9721E-01 
5, 1162E-01 
5.2505E-01 
5.3743E-01 
5.4874E-01 
5. 5895E-01 
5, 6803E-01 
5.7597E-01 
5.8275E-01 
5.8837E-01 
5.9281E-01 
5.9608E-01 
5. 9819E-01 
5.9913E-01 
5, 9892E-01 
5, 9758E-01 
5.9513E-01 
5.9158E-01 
5.8696E-01 
5, 8130E-01 
5.7463E-01 
5.6698E-01 
5.5839E-01 
5.4890E-01 



SF 

0, u 

•1, 9186E-04 
-1.6383E-03 
-4. 9661E-03 

-2.6190E-02 
-3.6636E-02 
*4. 8559E-02 
-6. 1829E-02 
-7.6317E-02 
-9. 1894E-02 
-1.0843E-01 
-1.2581E-01 
*1. 4389E-01 
-1.6256E-01 
•1. 8 170E-01 
*2, 0 1 19E-01 
-2.2091E-01 
-2, 4076E-01 
*2, 6063E-01 
*2. 8042E-01 
*3. 0002E-01 
-3.1935E-01 
*3, 3832E-01 
-3, 5684E-01 
-3.7483E-01 
-3. 9222E-01 
-4. 0895E-01 
-4. 2494E-01 
-4, 4014E-01 
-4. 5449E-01 
*4. 6795E-01 
-4, 8047E-01 
-4.9202E-01 
•5.0256E-01 
-5. 1206E-01 
-5.2051E-01 
• 5. 2788E-01 
-5.3415E-01 
•5, 3933E-01 
-5. 4339E-01 
•5.4635E-01 
*5, 4821E-01 
-5.4896E-01 
•5.4863E-01 
•5.4722E-01 
-5.4476E-01 
■5, 4127E-01 
*5. 3675E-01 
*5, 3126E-01 
•5. 2480E-01 
*5. 1742E-01 
-5.0915E-01 
•5. 0001E-01 



s *h 

0,0 

-4.2725E-06 

-7.8408E-05 

-3.7232E-04 

-kSsIii-BI 

-4. 1039E-03 
-6.6987E-03 
-1.0113E-02 
-1.4404E-02 
-1.9608E-02 
-2. 5750E-02 
-3.2838E-02 
-4.0866E-02 
-4.9819E-02 
-5.9671E-02 
-7.0385E-02 
-8. 1918E-02 
-9.4222E-02 
-1.0724E-01 
-1.2091E-01 
-1.3516E-01 
-1.4993E-01 
-1.6516E-01 
-1.8075E-01 
-1.9664E-01 
-2.1276E-01 
-2.2903E-01 
-2.4538E-01 
-2.6172E-01 
-2, 7800E-01 
-2.9413E-01 
-3.1005E-01 
-3.2569E-01 
-3.4098E-01 
-3, 5587E-01 
-3.7030E-01 
-3 ,842 IE-01 
-3.9754E-01 
-4, 1026E-01 
-4.2231E-01 
-4.3365E-01 
-4.4425E-01 
-4.5406E-01 
-4.6307E-01 
-4.7123E-01 
-4.7854E-01 
-4.8496E-01 
-4.9049E-01 
-4.9511E-01 
-4.9882E-01 
-5.0160E-01 
-5.0346E-01 
-5.0439E-01 
-5.0441E-01 



Table 3* lb. Computer output data for C, , Sp, and S^g. 
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00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

) . 



•0328k 00 
•0430E 00 
•0526E 
.06156 



Sk 3 



00 

00 



.06986 00 
• 0775E 00 
•08456 00 

• 091 0E 00 
•0969E 00 
.10216 00 
• 10696 00 
• 111 IE 00 
• 11486 00 

• 1 179E 00 
•1206E 00 

00 
00 
00 



• 1228E 
• 12466 
•1259E 



• 1268E 00 
• 1273E 00 
• 12746 00 
• 1271E 00 
' 1666 00 
>566 00 
144E 00 
•1229E 00 
00 
00 

•1168E 00 
00 
00 
00 
00 
00 
00 
00 
00 
00 



:tntr E 



• 1143E 
•1116E 
• 1087E 
.10566 
.10246 
• 09916 
•0956E 
•09206 
• 0883E 
•0845E 00 
•0806E 00 
•0767E 00 
•0728E 00 
• 0688E 00 
•0647E 00 
•0607E 00 
• 0567E 00 
•0527E 
• 0486E 
•0446E 
•0407E 00 
• 0368E 00 
•0329E 00 
•0291 E 00 
•0253E 00 
•0216E 00 



5.3854E-01 
5.27366-01 
5.1539E-01 
5.02696-01 
4.8929E-01 
4.75256-01 
4.60606-01 
4.4540E-01 
4* 2969E —01 
4, 13526-01 
3.96946-01 
3.7999E-01 
3.6273E-01 
3.4520E-01 
3.27446-01 
3.09516-01 
2.91446-01 
2.73296-01 
2.5509E-01 
2.36906-01 
2.18746-01 
2.00666-01 
1.82706-0 



=8 = 



00 

00 

00 



1.64906 
1.47286 
1.2990E-01 
1 • 12 776-01 
9. 59406-02 
7.94266-02 
6* 32606-02 
4.7470E-02 
3.2081E-02 
1.71146-02 
2.59216-03 
-1.14656-02 
-2.50406-02 
-3.8116E-02 
-5.06786-02 
-6*2716Er-02 
-7.42166-02 
-8.5169E-02 
-9.55676-02 
-1.0540E-01 
-1.14686-01 
-1.2338E-01 



39076-01 
-01 



-1.46066- 



-J.3151E-01 

-1* 52476-01 
-1.58336-01 
-1. 63626-01 
-1.6836E-01 
-1.7255E-01 
-1.7620E-01 
-1.7932E-01 



Sf 

4. 90066-0* 
4. 79326-0 
4.67846-0 
4.5566E-0 
4.42826-0 
4* 2937E— 0! 
4* 15356-0 
4.00796-0 
3.85766-0 
3.70296-0 
3.5443E-0 
3.38236-0 
3.2173E-0 
3,04976-0 
2. 8800E-0 
2.7086E-0 
2.53606-0 
2.36256-0 
2.18876-0 
2.01486-0 
1.84146-0 
1.66876-0 
1. 49726-0 
1* 32726-0 . 
1. 15906-0 
9.93026-02 
-8.29526-02 
-6.68826-02 
-5.11216-02 
•3* 56916-02 
-2.06226-02 
-5.93576-03 
8.3437E-03 
2.21996-02 
3. 56096-02 
4.85546-02 
6.10266-02 
7.30066-02 
8* 4481c— 02 
9. 5444E-02 
1. 0588E-01 
1. 1579E-01 
1.2516E-01 
1.3399E-01 
1.4227E-01 
1. 5001E-01 
1.57206-01 
1*63856-01 
1* 6995E-01 
1* 7551E-01 
1.8053E-01 
1. 8502E-01 
1.8899E-01 
1* 9245E-01 
1.95396-01 



Sth 

-5.0352E-01 

-5.01 72E-01 

-4.9904E-01 

-4.95496 -01 

-4.9108E-01 

-4.8584E-01 

-4.79806-01 

-4.7296E-01 

"4*65376-01 

-4* 57056-0 

■4*48036-0 

-4* 3834E-0 

-4.2802E-0 

■4*17106-0 

-4*0561E-0 

-3*93596-01 

-3.8107E-01 

-3.68106-01 

-3*54716-01 

-3.4094E-01 

-3.26836-01 

-3.12426-01 

-2.97746-01 

-2* 8283E-01 

-2.67736-01 

-2.5247E-01 

-2.3710E-01 

-2.21656-01 

-2.0616E-01 

-1.90666-01 

-1.75186-01 

-1.5976E-01 



■1.44426-01 
-1.29216-01 
1.14156-01 
•9. 9269E-02 
-8.4595E-02 
-7.01536-02 
-5.59716-02 
-4.2070E-02 
-2.84746-02 
-1.52056-02 
-2.28156-03 
1.02786-02 
2.24566-02 
3.42376-02 
4.56076-02 
5.65526-02 
.70626-02 
.71266-02 
8.67366-02 
9.58836-02 
1.04566-01 
1.12776-01 
1.2050E-01 



♦ 



mtinued ) 



40 



IV. EXAMPLES OF DIFFERENTIAL EQUATION SENSITIVITY ANALYSIS 



A. INTRODUCTION 

A general method of determining the sensitivity functions 
for the coefficients of a second-order differential equation 
is derived for a mass, spring, and damper system. The 
procedure is then applied to the study of a sixth-order 
differential equation, the solution of which is initially 
unstable. The coefficients are then adjusted to yield a 
stable solution using the results of sensitivity analysis. 

B'. SECOND-ORDER DIFFERENTIAL EQUATION SENSITIVITY TO 
COEFFICIENT VARIATIONS 

The elementary mass-spring-damper system shown in Fig. 

4.1 is described by the following differential equation 

MX + BX + KX = u(t) (4.1) 

Taking the Laplace transform and solving for X(s) gives 

X (s) = — (4.2) 

Ms +Bs+K 

For convenience, let M = 1. Applying (2.1) to (4.2) yields 

S K = U(s) - - 2 = X ( s) (— ^ ) (4.3) 

(s +Bs+K) s +Bs+K 

By comparing (4.3) to (4.2) it can be seen that the sensi- 
tivity model is a replica of the original system with the 
controlled output as the input to the model. The sensitivity 
coefficient is measured at the point corresponding to the 
feedback signal. 
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Fig. 4.1. Block diagram of mass- spring-damper system. 




-S 



K 



-Sb 



Fig. 4. 2. Second-order differential equation sensitivity 
model . 
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Similarly , 



S = - Bs = X ( s) ( — ) (4.4) 

(s 2 +Bs+K ) 2 S 2 +Bs+K 

Figure 4.2 is a block diagram of the complete sensitivity 
model which shows clearly that and Sg are comparable to 
the feedback gain sensitivities as derived in Section II 
and III. 

The solution of (4.1) for a unit-impulse input is 
1 _ B 

x(t) = e 2 sin \T~^ B 2 t (4.5) 

/ K S 2 K ' r- 

K T~ 



Figure 4.4 is the graph of x(t) and the sensitivity coeffi- 
cients S K and Sg for B = 0.6 and K = 0.6. The plot shows 
that the system is more sensitive to changes in K than B 
as verified by (4.5) . Table 4.1 is the computer output data 
for X/ S R , and Sg . 

C. SIXTH-ORDER DIFFERENTIAL EQUATION SENSITIVITY TO COEFFI- 
CIENT VARIATIONS 



Consider the differential equation 



d 6 x ,d 5 x d 4 x d 3 x d 2 x dx 

a + b + c + d + e + f — + g x 



(4.6) 

u(t) 



dt f 



dt 



dt 



dt ' 



dt ‘ 



dt 



where nominally a = 1, b = 3.4, c = 11.4, d = 9.4, e = 7.4, 
f = 0.575, g = 0.452. The solution when u(t) is a unit-step 
can be determined by DSL using only one block. For the 
sensitivity model the equation must be integrated in six 
steps to isolate eaph sensitivity coefficient. The procedure 
used for the second-degree equation can be expanded to give 
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X(8) 
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Fig. A. 3. Sixth-order differential equation sensitivity model. 




TIME (SEC ) 

Fig. 4. 4. Second-order system response and sensitivity 
functions . 
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Fig. A. 5 . Sixth- order unstable system response and 
sensitivity functions. 

(a) . Response and sensitivity functions for 

d and f. 

(b) . Sensitivity functions for b and o. 

(c) . Sensitivity functions for e and g. 




Fig. A. 5(a). X, Sfl , and Sf. 



46 



2 



1 



0 



-1 



-2 













_ JL 


^ s ° 

r\ 










3c: — V — 




?GJ 
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TIME ( SEC ) 



Fig.A.5(b). St, and S c . 




Fig.4.5(c). S e and Sg. 
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Fig. 4. 6. Sixth-order stable system response and 
sensitivity functions. 

(a) . Response and sensitivity functions for 

d and f. 

(b) . Sensitivity functions for b and c. 

(c) . Sensitivity functions for e and g. 




TIME (SEC) 



Fig. 4. 6 (a). X, S^ , and Sf. 
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1.0 



o.5 




-0.5 



-l.o 

Fig. A. 

4 

2 

0 

-2 



• TIME (SEC ) 

Fig.4 .6 (c ) . S e and S g . 



• TIME (SEC ) 
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a model similar to Fig. 4.2. The expanded sensitivity model 
is shown in Fig. 4.3. The solution of (4.6) is unstable as 
can be seen in Fig. 4.5(a) . The sensitivity functions 
graphed in Fig. 4.5(b) indicate that the inner-loop feed- 
back gain parameters b and c influence the output relatively 
little. Going from the coefficient of the highest order 
term to the lowest, the influence increases but it is ob- 
tained with an increasing time delay. 

The problem of stabilizing the solution may be solved 
using sensitivity analysis by noting that is out-of -phase 
with X. By increasing the coefficient f, the response 
becomes stable as shown in Fig. 4.6(a) for f = 1.8. Any 
further adjustment of the performance can be achieved by 
referring to the sensitivity functions for the stabilized 
system as shown in Fig. 4.6. 
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TIME 

0.0 

1 . 0000^-01 

2.0000E-01 

3.0000E-01 

4 • OOOOE-0 1 

5 .OOCOE-O 1 
6.0000E-01 
7. 0000 E- 01 
8. 0000 E- 01 
9. 0000 E- 01 
1.0000E 00 
1.1000E 00 
1.2OC0E 00 
1.3000E 00 
1.400CE OC 
1.5000E 00 
1.60C0E OC 
1.70C0E 00 
l.BOGOE 00 
1.9000F OC 
2.0000E 00 
2 . 1000 p OC 
2.2000E 00 
2.3000E 00 
2.4000E 00 
2.5000E 00 
2.6000E OC 
2.7OC0E 00 
2.8000E 00 
2.9000E 00 
3.0000E 00 
3.1000E 00 
3.2000E OC 
3.3000E OC 
3.4000E 00 
3.5000E OC 
3.6000E 00 
3.7000E 00 
3.8000E OC 
3.9000E 00 
4.0000E 00 
4.1000E 00 
4 .2OC0E 00 
4.3000E 00 
4.4000E OC 
4.8000E 00 
4.60COE 00 
4.70COE 00 
4 . 80 COE 00 
4.90COE 00 
8 ,0000 c 00 
5.1OC0E 00 
5.2OC0E 00 
5 • 30 COE 00 
5.4OC0E 00 



X 

C. 0 

A . 8990E- 03 
1 • 91 85E-02 
A . 22 28E-02 
7.3336E-C2 
1.12Q0E-C1 
1. 5742E-01 
2 . 0888E-01 
2.6601E-C1 
3 .27086-01 
3.939AE-C1 
A . 63 56E-C1 
5. 3613E-C1 
6 • 1 106E-0 1 
6.8776E-C1 
7 . 65 70E-C1 
8 .44 3AE-C1 
8 . 2319E-C1 
1.0018E OC 
1.0797E OC 
1.1564E 00 
1.2317F 00 
1. 305 IE 00 
1.3763E 00 
1.4451E CO 
1.5112E 00 
1.5743E CC 
1.63A3E 00 
1.6909F OC 
1.7A42E OC 
1.7939F CC 
1.8399E OC 
1.8823E 00 
1.92C9E CC 
1.9558E 00 
1.9870E 00 
2.01A5E CC 
2.038AE CC 
2.C587E CC 
2.0756E CO 
2.0891E CC 
2.0984E 00 
2.1065E CC 
2.11C7E CC 
2.1120E CO 
2.11C7E CC 
2.1068E OC 
2.1006E 00 
2.09235 00 
2.082CE CC 
2.06<?8E 00 
2.0560E 00 
2.0A08E OC 
2.0242E 00 
2.0066E 00 



SENSB 

- 0.0 

-9. 6995E-05 
-7.5187E-C4 
-2.4563E-03 
-5.63C1P-C3 
-1 • 0622E-02 
-1.7713E-C2 
-2. 71 1 AE-02 
-3. 8976P-C2 
-5.3385E-02 
-7. 0372E-C2 
-8.9912E-02 
-1.1193F-C1 
-1 . 3632E-C1 
-1.6290«=-01 
-1 • 9150E-C1 
-2. 21 87E-C1 
-2. 5377E-01 
-2. 86 Q 2E-01 
-3.2102E-01 
-3.5577P-C1 
-3. 9085E-C1 
-A. 2 594E-01 
-A.6072E-C1 
-4. 9487E-01 
-5.2807E-01 

- 5. 6 0C3F-01 
-5. 9C45E-01 
-6. 19C5E-01 
-6.4558E-C1 
-6.698 IE-01 
-6. 9151E-01 
-7.1049E-01 
- 7 * 2659F-01 

- 7. 3966E-01 
-7.4858E-C1 
-7. 5628E-C1 
-7. 5968E-01 
-7. 5 974E-01 
-7.5646E-01 
-7.4985E-01 
-7. 3995E-01 
-7.2682E-01 
-7. 1055E-C1 
-6.9125E-01 
-6. 6806E-C1 
-6.441CE-01 
-6. 1656E-01 
-5.8661E-01 
-5. 5445F-01 
-5. 2027E-01 
-4. 843CE-01 
-4.4676F-01 
-4. 0788E-01 
-3.6789E-01 



SENSK 

- 0.0 

-2.4401 E-06 
-3. 8079E-05 
-1 . 8789E-04 
-5.7840E-04 
-1. 3746E-03 
-2. 7729E-03 
-4. 9943E-03 
-8. 2778E-03 
-1.2874E-02 
-1.9041E-02 
-2.7034E-02 
-3.7106E-02 
-4. 95 00E-02 
-6.4443E-02 
-8.2148E-02 
-1. 0280E-01 
-1.2657F-01 
-1.5360E-01 
-1 • 8399E-01 
-2. 1783E-01 
-2.5516F-01 
-2. 9600E-01 
-3.4034E-01 
-3 . 8812E-01 
-4. 3928E-01 
-4.Q370E-01 
-5. 5123E-01 
-6. 1 173E-01 
-6.7493E-01 
-7.4077E-01 
-8.0885E-01 
-8. 7898E-01 
-9.5086E-01 
-1.0242E 00 
-1.0987E 00 
-1.1740E 00 
-1.2498F 00 
-1.3258E 00 
-1.4017E 00 
-1.4770E 00 
-1.5515E 00 
-1.6249E 00 
-1.6968E 00 
-1.7669E 00 
-1.8349E 00 
-1.9006E 00 
-1.9637E 00 
-2.0239E 00 
-2. 0809 E 00 
-2.1347E 00 
-2.1849E 00 
-2.2315E 00 
-2.2742E 00 
-2.3130E 00 



Table 4.1. Computer output data for X, Sg, and Sr* 
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TIMI 




X 




SEN 58 SENSK 


5 .5000*' 


OC 


1.998C5 


CC 


-3.27C3F-C1 -2 .34795 


00 


5 . 60 C0 C 


00 


1.5686E 


CC 


-2. 8555E-C1 -2.3784 c 


00 


5 .70005 


00 


1.8485E 


CC 


- 2.4366E-C1 -2.4 0495 


00 


5 .80005 


00 


1.828C5 


CC 


-2.0162E-C1 -2.4271E 


00 


5.90C05 


OC 


1. 90725 


CC 


-1.5965E-01 -2.44525 


00 


6.00005 


OC 


1.88625 


OC 


— 1.1799F— Cl -2.4591E 


00 


6.10005 


00 


1.86515 


CO 


-7.681 55-C2 -2.4688E 


00 


6.200CF 


OC 


1 . 944 IF 


CO 


-3.63e2F-C2 -2.47455 


00 


6.300CE 


OC 


1.8233F 


OC 


3. 12685-07 -2.4761E 


00 


6 .40005 


OC 


1.80285 


CC 


4.152CE-C2 -2.4739E 


00 


6.50005 


OC 


1.78275 


CC 


7.8616F-02 -2.4679F 


00 


6 .60005 


00 


1.76315 


CC 


1 • 1425 P -Ci -2.45825 


00 


6.70005 


OC 


1 . 744 IF 


CC 


1.48275-C1 -2.4451 E 


00 


6.80005 


00 


1.72575 


OC 


1 . 9C54F-C1 -2.4296E 


00 


6.90005 


OC 


1 . 709CE 


CC 


2. 1C92F-C1 -2.4090E 


00 


7 .00005 


OC 


1 .69125 


CC 


2.3931F-C1 -2.3865E 


00 


7.10C05 


00 


1.6752F 


CO 


2.656 1 5-C 1 -2.36125 


00 


7 .20005 


OC 


1.660CE 


CC 


2. 0C 74F-C1 -2.3334c 


00 


7.30005 


00 


1.6458F 


CC 


3.11635-01 -2.3033E 


00 


7.40CC5 


OC 


1.63265 


CC 


3. 3123E-C1 -2.27125 


00 


7.50005 


00 


1.6203F 


CO 


3.49515-C1 -2.23725 


00 


7.60CC5 


OC 


1 .60895 


CC 


3 . 6 344 F- Cl -2.2016E 


00 


7.70005 


00 


1.5986E 


CC 


3.76C15-C1 -2.1646E 


00 


7.80005 


00 


1.58935 


CC 


3 • 86 23 F -C1 -2.12645 


00 


7 .90005 


OC 


1.5910F 


CC 


3. 94125-C1 -2.0874E 


00 


8.0000E 


00 


1.5736F 


CO 


3.99725-C1 -2.0477E 


00 


8 . 10 005 


OC 


1.56735 


CC 


4.C3C6E-C1 -2.00755 


00 


8 .20005 


00 


1.56185 


OC 


4.C42C5-C1 -1.9671E 


00 


8 . 30CC C 


OC 


1.5574F 


CO 


4. 0 3 22 P - Cl -1.9268E 


00 


8.40005 


00 


1.55385 


OC 


4.0C17E-01 -1.8666E 


00 


8.50005 


OC 


1.551CE 


CC 


3.95165-01 -1.84685 


00 


8.60005 


00 


1.5481E 


CC 


3.8828F-C1 -1.80765 


00 


8.70C05 


00 


1.5480E 


CC 


3 . 79625-01 -1.76925 


00 


9.80CC5 


00 


1.54775 


CC 


3.6930E-C1 -1.7317E 


00 


8 ,900C C 


00 


1.54805 


CC 


3.57435-01 - 1 . 6 C 5AE 


00 


9.00C0F 


OC 


1.5481E 


CC 


3.441 35-C1 -1.66035 


00 


«.10005 


00 


1.55C7E 


CC 


3. 2952 5- Cl -1.62665 


00 


9.20CCE 


OC 


1.5529E 


OC 


3.1373E-C1 -1.59445 


00 


9.3000E 


00 


1.5557E 


CC 


2 . 96895— Cl -1.56305 


00 


9.40005 


00 


1.5590E 


CO 


2.7913E-C1 -1.5351E 


00 


9 .50005 


00 


1.56275 


CC 


2.6C575-C1 -1.50815 


00 


9.60005 


00 


1.5667E 


OC 


2.413 R 5-C1 -1.48305 


00 


°. 70005 


OC 


1.57125 


CC 


2. 216CF-C1 -1.4598E 


00 


8 . 80C0 C 


00 


1.5759E 


CC 


2.C144F-01 -1.43875 


00 


8.90CCE 


00 


1.5809E 


CC 


1.81 CC C -C1 -1.41055 


00 


1 .00005 


01 


1.5860E 


OC 


1.6C4C5-C1 -1.40255 


00 


1.01005 


01 


1. 5«14E 


CO 


1 . 39765-C1 -1.3875E 


00 


1 .0 2005 


01 


1.59695 


CO 


1.192CE-C1 -1.37455 


00 


1.03005 


01 


1.6024F 


CO 


5.8811F-02 -1.36365 


00 


1 .04005 


01 


1.60P1F 


CC 


7.87135-02 -1.3547E 


00 


1 .05005 


01 


1.6137F 


OC 


5 • 8999F-02 -1.34795 


00 


1.06005 


01 


1.6193E 


00 


3.97635-02 -1.3429E 


00 


1.07005 


01 


1.6249E 


OC 


2.1C90E-C2 -1.3399E 


00 


1.08005 


01 


1.6303E 


OC 


3.06125-03 -1.3387E 


00 


1 .09005 


01 


1.6357E 


CC 


-1.42525-02 -1.33925 


00 



Table 4.1. (Continued) 
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V. NONLINEAR SYSTEM SENSITIVITY 



A. INTRODUCTION 

While it is true that no physical systems are exactly 
linear, often adequate linear approximations can be used. 

In these cases, the Kokotovic [6] method of sensitivity 
analysis as developed in the preceding sections of this 
thesis can be applied. There are many situations, however, 
for which linear representations are invalid and nonlinear 
characteristics must be considered. 

Figure 5.1 is a block diagram of a nonlinear system 
containing a linear element represented by a transfer function 

K , where D = d/dt is the differential operator, and a 

nonlinear element, N, whose output is denoted by f(e) . D is 
used instead of s in the linear transfer function because the 
Laplace transform and its inverse are generally not suitable 
for nonlinear analysis. 

The sensitivity to gain-constant variations is developed 
for the general nonlinear system and then applied to a 
specific system containing a saturation element as the non- 
linear component. 

B. SENSITIVITY COEFFICIENT FOR A GAIN FACTOR 

In order to designate the dependence of the variables in 
Fig. 5.1 the following notation is used: 
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Fig.5-1* Block diagram of nonlinear system. 




Fig. 5. 2. Sensitivity model of nonlinear system. 
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r 




r ( t) 




e 




e ( t ,K) = r (t) - c(t,K) 




f (e) 


= 


output of nonlinear element N 




c 


= 


c(t,K) = r (t) - e ( t ,K) 




X 


= 


X ( t,K) 




P(D) 


= 


,n ,n-l ^ 

a d + a n a . + ••• + ax Jr + a o 

n dt n n_1 dt n_1 dt 


(5.1a) 


Q(D) 


= 


,n ,n-l * 

b -rr- + b , r- + • • • + bi 3-r- + 

n dt n-1 ^ n-1 dt 


: (5.1b) 



From the definition of gain-constant sensitivity coeffi- 
cient (2.1) 



3 c _ y 3 e 

K ” 3K/K “ 3K 



(5.2) 



Using the differential equation for the linear element of the 
system gives 



KP(D) f (e ( t,K) ) = Q (D) c(t,K) = Q(D) [r ( t) -e ( t ,K) ] (5.4) 



Taking the logarithmic derivative, 3/3 in K = K 3/3 K, of 



(5.4) and 
KP (D) 



substituting (5.2) yields 
f (e ( t , K ) ) + K 2 P (D) ^ ) ■ s 

_k Q( d)|#1^kI 



(5.5) 



Rearranging and simplifying (5.5) gives 



S 



K 



£e = KP (P) f /_» , K 2 P (D) 3f (e) 

3K Q(D) ' 1 Q (D) 3K 



(5.6) 



Figure 5.2 is a block diagram of the proposed sensitivity 
model for the nonlinear system. It will be shown that the 
describing equation of this model is identical to (5.6) . 
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From Fig. 5.2 



c _ 3c _ , 

S K 3K/K C_b 

Substituting (5.2) in (5.7) yields 



(5.7) 



(5.8) 



Applying the transfer function for the linear element of 
the model gives 



-K 2 P (D) a | - ( — ||- = Q (D) b = Q ( D) c + KQ (D) |f (5.9) 

Substituting (5.4) in (5.9) for Q(D)c and solving for S^ 

Ja 

as defined by (5.2) gives 



K 



K3e 

cTk 



= K 



P(D) 

qW 



f (e) + K 



2 P(P) 3f (e) 
Q (D) 3K 



(5.10) 



Equations (5.10) and (5.6) are the same and, hence, the model 
generates the gain-constant sensitivity coefficient of inter- 
est. The sensitivity model block represented by ^ e — is the 
instantaneous slope of the nonlinear function f (e) . In many 
nonlinear systems a piecewise-linear approximation can be 
utilized. The sensitivity model will be the same as that of 
a linear system except that the gain changes as a function of 
the error signal in the original nonlinear system. An example 
for a piece wise-linear saturable element will make this more 
meaningful . 



C. APPLICATION 

The nonlinear component used as an example will be a 
limiter with a saturation characteristic as shown in Fig. 5.3. 
The linear plant used in Section II is used for comparison 
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of the time responses and gain sensitivity coefficients. 

The transfer function in differential operator notation 
becomes 

KP (D) _ 50 

Q ^ D 3 + 10D 2 + 24D 

The gain-constant sensitivity model as shown in Fig. 5.2 
requires a piecewise-linear function described by Fig. 5.4. 
It can be seen that ^ — is the slope of the nonlinear 

function f(e) and may be handled by logic statements in the 
DSL computer program as given in Appendix D. 

Figure 5.5 shows the time response and the sensitivity 
function for a unit-step input with a limited amplitude of 
error signal set at 0.1. The nonlinear element is saturat- 
ed until 4.8 seconds. During this period of saturation the 
gain representing the nonlinear element in the sensitivity 
model is zero, and therefore, the sensitivity coefficient 
and the output of the nonlinear system are one and the same 
as shown in Fig. 5.5 for the first 4.8 seconds. Figure 5.6 
shows the response and sensitivity for the limited amplitude 
of pi = -1.0 and p 2 = 1.0. This lets the system perform in 
a completely linear fashion and the effect of the nonlinear 
element may be realized by comparison with Fig. 5.5. The 
nonlinear response is slowed down, the maximum sensitivity 
to gain is greater but delayed by 3.8 seconds; also, the 
peak overshoot is greatly reduced. 

To verify the sensitivity coefficients for gain vari- 
ations, check points were chosen in the saturation region, 
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f(e)=P!, 
f (e )=p 2 , 
f (e )=e , 



?i-S*5»3* Characteristics of nonlinear N. 




If (e ) = 1 
9 e 

2 f (e ) = 0 
d e 

d f (e ) - 0 
a e 



Fig. 5. A. Characteristics of sensitivity model 



e< Pi 
©>P2 

Pl^e < p 2 



, Pl<e < p 2 
, e <Pi 
, e > p 2 



nonlinearity. 
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1.5 



1.0 



0.5 



0.0 



-0.5 
Fig. 5. 

1.5 



1.0 



0.5 



0.0 



-0.5 
Fig. 5. 




5. Response and sensitivity for nonlinear system 




6. Response and sensitivity for unlimited linear 
system . 
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peak-positive, and peak-negative portions of the sensitivity 
curve as shown in Fig. 5.7 for pi = -0.2 and p z = 0.2. The 
approximate change in the response due to a fractional change 
of 0.1 in the gain parameter, 

AC - S„ AK/K 

was calculated and compared to the computer output for K = 55 
and K = 45 as given in the Table 5.2. The maximum error 
found was 0.5% for point 2 as presented in Table 5.1. 

Figure 5.8 shows the response for 14% variations in 
gain as was required for the linear system to constrain the 
peak overshoot to within + 0,05. The peak overshoot for 
the nonlinear system was constrained to + 0.014 for the same 
gain variations. This can be predicted from Fig. 5.8 since 

' C, *• 

at the time of peak overshoot, the sensitivity function is 
much nearer zero than the sensitivity function for the linear 
system shown in Fig. 5.7. 
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1.5 




J TIKE (SEC ) 

FiS. 5.7 • Response and sensitivity for nonlinear system 
with p]_=-0.2 and P2= 0.2. 



Pt. 


T 


c 50 


S K 


C+AC 


c 55 


0 

1 

> 

O 


C45 


1 


1.5 


.45215 


.45215 


.49737 


.49736 


.40694 


.40693 


2 


2.5 


.86763 


.84752 


.95238 


.94755 


.78288 


.78127 


3 


4.0 


1.0177 


-.21111 


.99659 


.99673 


1.0388 


1.0351 



Table 5.1. Comparison of actual outputs and outputs 
predicted from sensitivity coefficients. 
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8. Gain tolerance of 14$ for overshoot of ±0.014. 
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K=50 



T IMF 




SENSK 


n 17 


0.0 




c.o 


0.0 


l.OOOOE- 


01 


1. 304 1E-C3 


1.3C416 


2.00009- 


01 


8. 30249-C3 


8.3024F 


3.00009- 


01 


2.25529-02 


2.2552F 


4.00C0F- 


01 


4.3547E-02 


4.35479 


5.00006- 


01 


7.0099E-02 


7.0C99E 


6.00006- 


01 


1. 00949-01 


1.0C94F 


7.00009- 


01 


1.3498E-C1 


1.3498E 


8.0000E- 


01 


1 .71326-01 


1.71326 


9.00006- 


01 


2 • 0930E- Cl 


2. 093C C 


l.OOOOE 


00 


2.4843E-C1 


2.4843F 


1.10006 


00 


2. 8837E-C1 


2.8837E 


1.2000E 


OC 


3.2886E-C1 


3.2886F 


1.3000E 


OC 


3.6972E-C1 


3.6972 c 


1 .4000 c 


00 


4. 1085E-C1 


4.1085F 


1.5000E 


00 


4.521 5E-01 


4.5215E 


1.6000E 


OC 


4.9357E-C1 


4.9357F 


1 .70009 


00 


5 . 35079-C1 


5.35C7F 


1.8000E 


00 


5 .76626-01 


5.76629 


1 .9000 c 


OC 


6 • 1821E-C1 


6. 1 82 IF 


2.0000E 


00 


6 . 5983E- Cl 


6.5983F 


2.1000E 


00 


7.0146E-C1 


7.0146F 


2 .20009 


00 


7.431 1E-C1 


7.431 IF 


2. 3000E 


CC 


7 . 84769-C1 


7. 84769 


2.40006 


OC 


8. 24946-01 


8.26419 


2.500CF 


CC 


8.4752E-C1 


8.6763E 


2.6000E 


CO 


8 . 3906E-01 


9.072C«= 


2 .7000 c 


OC 


7.9771P-C1 


5.4366F 


2.80009 


OC 


7 . 2738E-C1 


9.75779 


2.9000E 


OC 


6.34 57E-C1 


1.CC27F 


3.00009 


00 


5 .26696-01 


1.C2399 


3.100CE 


00 


4.11136-01 


1.0395E 


3.2000E 


00 


2.94686-01 


1.0496F 


3.30009 


OC 


1 . 8319F-01 


1.C54PF 


3.40009 


CC 


°. 13446-02 


1. 05 c 89 


3.500CF 


CC 


- 7 . 3481E-C3 


1.C5339 


3.60009 


00 


-8.061 69-C2 


1.C4839 


3.70009 


00 


- 1 . 37326- Cl 


1.C416F 


3.8000P 


00 


- 1.77369-C1 


1.C33PF 


3.9000E 


00 


-2 . 0147E-C1 


1.02569 


4.00C09 


00 


- 2 . 1 1 1 IE- Cl 


1.0177F 


4.1000E 


OC 


-2.0823E-C1 


1.C1C3E 


4.20009 


00 


-1.95C9E-C1 


1.CC38F 


4.30009 


00 


- 1 .74126-01 


9.98459 


4.4000F 


OC 


-1.4774E-C1 


9.9425E 


4.5000E 


00 


-1.18 266-C1 


9.91249 


4 .60 COE 


00 


-8.77356-02 


9.8934F 


4.7000E 


00 


-5.79539-02 


9.8844F 


4 .80 009 


00 


-3 . 0345E-C2 


9.88399 


4.90009 


00 


- 5 . 9862E-03 


9.89C1F 


5 .OOCOF 


00 


1.4402E-02 


9.90146 


5.10C0E 


OC 


3.0435E-C2 


9.91619 


5.2OC0E 


00 


4 . 2025E-C2 


9.9327E 


5.3000E 


00 


4 • 93336-C2 


9.94979 


5.40C0E 


00 


5.27236-02 


9.9662F 



OIT OIT 

C.C O.C 

C3 1.4345F-C3 1.1737F-C3 

03 9.1326F-03 7.4721E-C3 

C2 7.4AC7P-C2 2.C297E-C2 

C2 4.79029-02 3.91936-C2 

02 7.71 C9F-02 6.3089F-02 

Cl 1.11C4E-C1 9.084PE-02 

01 1.48476-01 1 .21489-01 

01 1.8845F-01 1.541 9E-01 

Cl 2.3023F-C1 1. 8837F-01 

01 2.7328F-C1 2.2359E-C1 

Cl 3. 1 7 2 16-01 2.5953«=-Cl 

01 3.6174E-01 2.9597E-01 

Cl 4 • 067CF-C1 3.3275E-C1 

01 4. 5 1 93E-01 3. 6976E-C1 

Cl 4 . 9736E-01 4.0693^-01 

Cl 5.4292F-01 4.44219-C1 

01 5 • 8857F-C1 4. 91 56F-C1 

01 6.3428F-01 6.1896F-01 

Cl 6.80C3E-C1 5.56369-01 

01 7.2581F-C1 5.9385F-C1 

Cl 7.716 1E-C1 6.3132F-C1 

Cl 8. 1742F-C1 6. 688CE-01 

Cl 8.6294E-C1 7.0628F-C1 

Cl 9.C691E-C1 7.4377F-C1 

Cl 9.4755F-01 7.8127E-01 

Cl 9. 833CE-C1 8.1 877F-C1 

01 1.013CF CO 8. 56CCE-C1 

Cl 1.0361E CO 8. 92C5E-C1 

CO 1.0523F CC 9.2573F-C1 

CO 1.0621F CO 9. 5599E-C1 

CC 1.C461F CO 9. 82 C7F-01 

CO 1.0652F CO 1.0035F OC 

CC 1.06C3F CC 1.02C3E 00 

CO 1.0526E CO 1.C324E CO 

CO 1.0431F CC 1.04C3F CC 

CC 1.03289 CO 1.04446 CO 

CO 1.0224E CC 1.C453E CC 

00 1.C126F 00 1.0436E CO 

CO 1 . 004 C p 00 1.04CIE CO 

CO 9. 9673F-C1 1.C351E CO 

OC 9.911 2F-01 1.C294F 00 

00 9.9717F-C1 1.C234P 00 

Cl 9.8479E-C1 1.C174E CO 

Cl! 9.8382E-C1 1.0118F 00 

Cl 5.84C6F-C1 1.0067E CC 

Cl 9. 8526E-01 1.0024E CO 

01 9. 8714F-01 9.9881F-C1 

Cl 9.9947E-01 9. 96C7E-C1 

01 9. 92CCF-C1 9.9412E-01 

01 9. 9454F-C1 9.9289E-C1 

Cl 9.9693E-01 9.9229E-01 

Cl 9. 99C5F-C1 9. 9222F-01 

01 1.00C9F CO 9.92579— Cl 

01 1.0022F 00 9.9323F-C1 



Table 5.2. Computer output data for nonlinear system 
outputs and sensitivity. 
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T I M F 

5.5000E 00 
5.60CCE 00 
5 .70 COE 00 
5.80CCE 00 
8.90C0E 00 
6.0000E 00 
6.100CE 00 
6.200CE 00 
6.3000F 00 
6.4000E 00 
6.5000E 00 
6.600CE 00 
6.7000E 00 
6.6000F 00 
6.9000E 00 
7.0000E 00 
7.1000E 00 
7 . 20 COE OC 
7.3000E OC 
7.4000E 00 
7.5000E 00 
7.6000E 00 
7.7000E 00 
7 . flOOOE 00 
7 . 90 COE 00 
8.00COE 00 
8.1000E OC 
9.20C0P 00 
8.3000E 00 
8.400CE 00 
8.5000F 00 
8 .60 COE 00 
8.70C0E 00 
8 • 80 COE 00 
8.9000E 00 
9.0000E 00 
9.1000E 00 
9.20C0E OC 
9.3000E OC 
9.40C0F 00 
9.5000E 00 
9.60C0E 00 
9.7000E 00 
9.8000F 00 
9 . 90 COE 00 
l.OOOOE 01 
1.0100E 01 
1.0200E 01 
I.0300E 01 
1.0400F 01 
1.0500E 01 
1.0600* 01 
1.0700E 01 
1.0800F 01 
1.0900E 01 



Table 5.2. 



SENSK 

5.27C3E-C2 
4 . 9877E-C2 
4.4889F-C2 
3.9392E-02 
3 • 1001F-C2 
2.3272E-C2 
1 . 5682E-C2 
8 . 6 14 IE- 03 
2.3564E-03 
-2.8982E-C3 
- 7 . 044 8E-03 
-1.0060E-02 
- 1 • 19P4E-02 
-1.2908E-02 
-1.2962E-C2 
-1.2303E-C2 
- 1 . 1 0 94E-C2 
-9.5025E-03 
-7.6836E-03 
-5.7786E-03 
-3 • 9067E-03 
-2.1646E-C3 
-6.2460E-04 
6.6608E-C* 
1.6829E-C3 
2 .4202E-C3 
2 • 8889*- 03 
3.11 24E-C3 
3. 1239E-03 
2.961 5E- 03 
2.6658E-C3 
2.2785E-C3 
1.8387E-C3 
1 . 3797E-C3 
9.3007E-04 
5 • 1326E-04 
1.4567E-C4 
-1.6022E-04 
-4.0054E-C4 
-8.741 IE- 04 
-6.8283E-04 
-7.3338E-C4 
-7.3433E-C4 
-6.9433E-04 
-6.2269E-04 
-5.2989E-C4 
-4.2456F-C4 
-3.16A4E-C4 
-2.1058E-C4 
-1.1283E-C4 
-2.7299E-05 
4 . 33 33E-05 
9.8586E-C5 
1.3816E-C4 
1.628AE-C4 



K=50 

OLT 

9.9813F-01 
q.CQ44F-01 
1.00C5E CC 
1.0013F 00 
1.CC19F 00 
1.0023E CO 
1.0024F CC 
1.0024* 00 
1.0023F CC 
1.002CE 00 
1.CC17E OC 
1.0013E 00 
1.0C1CF OC 
1.00C6F 00 
1.00C3E OC 
1.0001F OC 
9. 9985F-C1 
9.996 9 C -01 
9.998 8F— Cl 
9 . 9 95 IE- Cl 
9.9949E-01 
9. 995CF-C1 
9 • 995 3F-C1 
9. 99*9F— 01' 
9. 9966F-01 
9. 9973F- 01 
9. 9981F-C1 
9. 99R8F-C1 
9. 9994F-C1 
9.9999E-01 
l.OOCCE 00 
1.0001F CC 
1.0CC1F 00 
1.0001E CC 
1.0001F 00 
1.C0C1* CC 
1.CC01E 00 
1.0001E CO 
1.0CC1F CO 
1.00C1F CC 
l.OCOCE 00 
l.OOOOE CC 
l.OOOOF 00 
l.OCOCF 00 
9. 9999F-C1 

9. 999 8F-C1 

9. 999 8F-C1 
9.9998P-01 
9.9C98F-01 
9. 9998F-01 
9.999PF-01 

9 . 999 8E-01 
9.9998E-C1 
9. 9999E-01 
9. 9999F-C1 



CC 

00 

00 

CO 

CO 

00 

CO 

00 

00 

OC 

CC 



K=55 

mi 
. 003 IF 
.0037* 

. C04 CF 
. 0C39F 
. 0036F 
. C031 F 
. C026F 
.0C1 C * 

. 001 3E 
. 0CC7E 
. 0002E 
9. 998CF-C1 
9.9946F-01 
9.9923E-01 
9.99C9E-C1 
9 • 99 03F-01 
9. 9904F-C1 
9 . 991 2E-01 
9.9923F-01 
9.9937--C1 
9. 9952F-C1 
9.9967E-01 
9.9982F-C1 
9.9994F-C1 
l.OOOOE 00 
l.OCCl* 
1.0CC2* 
1.00C2E 
1.00C2F 
1.CCC2E 
1.0CC2E 
1.C002F 
1.00C2E 
1. 00C1F 
1 . 0001 * 

1. COCCE 
l.CGOC* 
9.9999F-01 
9. 9997E-01 
9 • 99 95F-C1 
9.9994E-C1 
9 .9994F-01 
9 . 9994F— 01 
9. 9994F-C1 
9. 9995F-01 
9. 9996E-C1 
9.9997E-C1 
9.9998F-01 
9.9999F-01 
9.99 c<;e-c 1 
l.OOOC* 00 
l.OOOOE 
l.OOOOE 
l.OOOCF 
l.OOOCE 



CO 

CO 

CO 

00 

00 

CO 

00 

CO 

00 

00 

CO 

CO 



00 

CO 

CO 

CO 



K=45 

Oil 

9.941 1F-C1 
9.951 1F-01 
9.9616P-C1 
9.971 8E-C1 
9.9813E-C1 
9.9899E-01 
9. 9971F-C1 
1.00C3F CO 
1.00C7E 00 
1.0011F CO 
1.0012F 00 
1.0013F 00 
1.0013F 00 
1.CC13* 00 
1.C011F 00 
1 . 001 CF OC 
1.00C8F 00 
1.00C6E CO 
1.00C5F 00 
1.00C3F CO 
1.C0C1F 00 
l.OOCCE 00 
9. 9993F-01 
9. 998fcE-Cl 
9.99P1E-01 
9.9978F-01 
9. 9977 F- Cl 
9.9977E-01 
9. 9978E-01 
9. 898CF-C1 
9.9983E-C1 
9 . 9986F-C1 
9.9989F-C1 
9.9992F-01 
9. 9995E-01 
9. 9997F-01 
9 . 9999F-01 
.COOCE 00 
. OOOCE CO 
. OCCCF CO 
.000C C CO 
• OCCCE 00 
. OOCCF CO 
. CCOCE CO 
. OOOCF 00 
• OCCOF 00 
.OOOCF CC 
. CCCCE 00 
.OOOCF 00 
• CCOCF CO 
.OOOCF CO 
.OOOCE 00 
9. 9999E-C1 
9.9999F-01 
9. 9999E-01 



( Continued ) 
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VI. CONCLUSIONS 



The technique of generating sensitivity functions for 
a basic linear feedback control system using the concept 
of a sensitivity model is developed. This model* is a replica 
of the original system and is cascaded with it. All sensi- 
tivity functions are generated simultaneously with the time 
response by use of the IBM/360 Digital Simulation Program. 

Sensitivity of a differential equation with respect to 
its coefficients is investigated. Computer graphs show the 
differences in the influence of the coefficients on the 
solution of the differential equation. The coefficients of 
the lower-order terms have the greatest influence but the 
effect is obtained with the greatest time delay. 

A theoretical sensitivity model is developed for a non- 
linear control system. Instead of a replica of the original 
control system as required for the linear model, the non- 
linear element is replaced by a component with the instanta- 
neous gain equal the slope of the nonlinear function., An 
example illustrates the procedure for a nonlinearity due to 
a saturation element. 

Sensitivity analysis indicates that proper interpretation 
of the sensitivity function results in accurate estimates for 
changes in the system response due to parameter variations. 
Tolerance limits of parameter variations are determined and 
performance improved by application of sensitivity coefficients. 
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APPENDIX A 



DSL/360 DIGITAL SIMULATION PROGRAM FOR BASIC SYSTEM 



TITLE BASIC SEPVC SFNSITIVITY COE c FI C I ENT S 
INTEG RKS 

STORAG IC(3» ,NLIM(1),DEN<4) 

TABLE ICC l-3)=4*0 t MJM( I) =1 .0, DEN (1-4) =1 . , 10. , 2^. 1 0 . 
CONTRL FINTIM=10. ,DELT=.l 
PARAM IN = 1. ,F = I . ,K=50. 

* SOLUTION C«= BASIC SYSTEM FCR TIME RESPONSE 

ERRORl=IN-F*CUT 

OUT=TRKFR (0»3,IC ,NUM,DEN,K*ERR0R1) 

* SOLUTION CF MODEL BOR SENSITIVITY COEFFICIENTS 

* SENSK=SENSITI VITY TO FORWARD GAIN K 

* SENSF=SEN$ITIVITY TO FEEDBACK GAIN R 

SFNSK=CLT+$FNSF 

X = TRNFR(0»3»IC» NUM »DEN»K*S C NSK) 

SEN SF = X* ( -F ) 

PRINT .l.OUT.SFNSK ,SFNSF 

PREPAR .1 ,CL'T t SENSK,SENSF 

GRAPH SAME » T I MF ,CUT , SENS K , SENSE 

PRPLOT ONLY 

END 

STOP 
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APPENDIX P 



DSL/36C DIGITAL SIMULATION PROGRAM FOR MULTILOOP SYSTEM 



TITLE MULTUCCP SERVO SENSITIVITY COEFFICIENTS 
INTEG RKS 

CONTRL FINTIM=1C.0»DELT=.01 

PARA* Kl=40. ,K2=1. »K3 = 1. ,K4 = ?. t FGAIN=l., IN=1 . t TAUI = 4., .. . 
TAU4=1. 

* SOLLTICN CF MLLTILCCP SYSTFV FOR TIME RESPONSE 

errcpi=in-fgain*cut 

D = EXCI TM-K4)*TAU4 
FEDl=D 

FED2=RFALPL( 0. ,TAU4,D) 

FEDR=FEC1-FEC2 

EPROR2=FRRCRl+FEDB 

EXCIT = RFALPL(0. ,TAU1 ,K1*frroR2) 

GEN=P c ALPL(0. ,2. ,K2*EXCIT) 

0UT = PFALPL(0. ,1 . » K3*GE N ) 

* SOLUTION CF VCDEl FOR SENSITIVITY COEFFICIENTS 

SFN$K3=CLT-FGAIN*Z 
0 X= X* ( -K4 ) * TAU4 
FED XI =CX 

FEDX?=PFALPL(0. ,TAU4,CX) 

?FNSK4=FEDX1-FEDX2 
$CNST4=REALPL(0. t TAU4,S c NSK4) 

SFNSK1=SENSK3+SENSK4 

P=-SFNSK] 

Q=RFALPL(0. tTAUl ,P) 

SEN ST 1 =P-Q 

X=RFALPL(D. ,4. »K1*SENSK1) 

Y=REALPL(C. »2. »K2*X) 

Z=RFALPL( 0. *1. ,K3*Y) 

PRINT . 1 ,OLT , SENSK3 tSENSKl ,SENST1 ,S E NSK4 , SENST4 

RREPAR .1 ,0UT,SENSK3 1 SFNSK1 .SFNS T I , S E NSK4 , SENST4 

GRAPH SAME » T I ME ,CUT . SE NSK3 , SENSK I , SENS T 1 » $ ENS K4 , SENST4 

GRAPH SA^E.TIME ,SENSK] t SENSTl 

PR PLOT ONLY 

END 

STOP 
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APPENDIX C 



DSL/36C DIGITAL SIMULATION PROGRAM FOR SIXTH ORDER 
DIFFERENTIAL FQLATIOM 



TITLE SIXTH ORDER DIFFERENTIAL EQUATION SENSITIVITY MODEL 
I N T EG RKS 

CONTRL FINTIM=50. ,DELT=.05 

PAR AM A = 1 . ,B = 3.4 ,C=11.4,D=9.4,E=7.4,F = .'575 ,G=.45 2 

* SOLUTION OF DE FOR TIME RESPONSE 

TN=STFP(C) 

XD6=(-B*XD5I +(-C*XD4)+(-D*XD3)+ (-E*XD2 )+(-F*XDl ) .. . 

+ ( -G+ X ) ♦ I N 

XD5 = IMGRL(0.,XD6) 

XD4 = I N TGP L ( 0. » XD5 ) 

XD3 = 1 N TGR L I 0 . ,XD4> 

XD2=INTGPL( 0. ,X03) 

XD1=INTGPL(0.,XD2» 

X=INTGRL( 0. , XD1 ) 

* SOLUTION OF MODEL FOR SENSITIVITY COEFFICIENTS 

SFNSG=-G*$X 
SENSF=-F*SXD1 
SENSF=-F*SXD2 
SENSD=-D*SXD3 
SFNSC --C* SXD4 
SEN SB^-B* S XD5 

SX06=SENSB+SENSC+SENSD+SFNSF+S6NSF+SENSG+X 
$ XDS = I NTGR L ( 0. , SXD6) 

SXD^=IMGRL( C. t SXD5 ) 

SXD3=INTGPL(C. »SXD4) 

SXD2=I NTGRL ( 0. ,$XD3) 

SXOI=INTGRL( C. t SXD2> 

SX=INTGRL(C. »SXD1) 

PRINT .1 tXtSEN SB, SEN SCtSENSD, SENSE, SENS F,SENSG 

PREPAR .1 ,X, SEN SB, SEN SC,SENSD, SENSE, SENS F,SENSG 

GRAPH SAME, TIME ,X, SENSB,SE NSC, SENSD, SENSE, SENS F,SENSG 

PRPLOT ONLY 

END 

STOP 
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APPENDIX D 



DSl/360 DIGITAL SIMULATION PROGRAM FOR NONLINEAR STSTFM 



TITLF NONLINFAR SYSTEM SENSITIVITY 

INTFG RKS 
STORAG IC( 3 ) «B ( 4 ) 

TABLE IC( 1-3)=3*G,B( 1-4) »1.,10. ,24. ,0. 

CONTRL FINTIM=12. ,CELT=.l 
PARAM K = 50. ,Pl=-.2 ,P2 = .2 

* SOLUTION CF NCM.INFAO SYSTEM FOR TIME RESPONSE 

IN=STFP ( C ) 

PRRCR = TN-CL'T 
Y = LI MI T ( PI ,P2, ERROR) 

0UT=DECM3,IC,B,K*Y) 

* SOLUTION OF w OCE L FOR SENSITIVITY COEFFICIENTS 

* SENSK=SENSI TIVITY TO FORWARD GAIN K 

NOSORT 

IFIFPRCR.GT. P2. OR. ERROR. LT.PDGC TO 1 

M=OL’T-N 

GO TO 2 

1 M=C 

2 CONTINUE 
N=DFQN(3»IC»B»K*M) 

SFNSK=CIT-N 

PRINT .1»FRRCR»CLT tSENSK 

PRFPAR * 1 »0L T , SEN SK 

GRAPH SAMEtTIMF »CUT,SENSK 

PRPLOT ONLY 

END 

STOP 
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